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The nine Peactor Statics Modules are designed to 
introduce students to the use of numerical methods and digital 
computers for calculation of neutron flux distributions in space and 
energy which are needed to calculate criticality, power distribution, 
and fuel burnup for both slow neutron and fast neutron fission 
reactors. The last module, BS-9r :_ncludes a separate program, WANDIC 
(one-dimensicnal diffusion code with three energy-group 
representation), for the calculations of neutron spectra and group 
constants for use in fast breeder reactor calculations. It can be 
used with the criticality programs in PS-2 and BS-8. The reactor is 
divided into several concentric annular regions and it is assumed 
that those regions can te represented fairly well by their average 
properties. The previous modules show how these properties are found. 
The set of nine modules is intended to supplement textbooks and other 
lecture material generally available to students in their course 
work. It is assumed that students are familiar with elementary 
nuclear structure, neutron-nuclei interactions, and introductory 
material cn fission chain reactors. (Author/SK) 
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Reactor Statics Module, RS-9 

Multlgroup Diffusion Program Using an Exponential Acceleration Technique 
9.1 Introduction 

A reactor physicist or a nuclear engineer designing a power reactor needs 
to know the number of neutrons present within the system. This information 
is necessary for the calculation of the power distribution in a reactor from 
which further vital information on other related physical quantities needed 
to assure a successful operation of the reactor can be derived. 

Previous modules discussed a broad area of problems which are encountered 
in reactor physics design. The basic input for calculations consists of 
engineering Sjpncif ications and nuclear data. The former provide quantities 
like power density, fuel composition and a geometric description of the 
reactor. This information varies from one design to another. On the other 
hand, the nuclear data contains detailed cross-section information. 

This input starts a rather long string of calculations before one can 
obtain results in a form which is suitable for engineering application. 

The detailed treatment of space-energy physical phenomena in the entire 
reactor represents a very difficult task. Therefore, the problem is usually 
divided into several different stages and the flow of calculations passes from 
the detailed analysis of the neutron balance In an elementary fuel pin cell 
to the reactor as a whole consisting of assemblies which are represented by 
their average properties. In this procedure, the spatial and energy effects 
are reduced in detail; the former through the "homogenizatlon", the latter 
through the group condensation. Of course, care must be taken in these simpli- 
fications. For example, neutron balance must be conserved. 
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The code WANDIC (one-dimensional diffusion code with a three energy-group 
representation) which is described in this module can be used to perform the 
global reactor calculations • The reactor is divided into several concentric 
annular regions and it is assumed that we can represent thor>: *.^^.io' s fairly 
well by their average properties. The previous modules shcv ^ . ge proper- 
ties are found. 
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Features of WMDIC 

9.2 General Description 

The program solves the three group diffusion equations 

DgVV - 2^*^ + = 0 for g = 1,2,3 (9.2.1) 

where the source term consisting of fission neutrons appears in the highest 
energy group only. The source terms in the lower groups are due to scattered 
neutrons from the adjacent upper group. The source S Is given as 

6 3 

= E (vl.)^*^ + ^^"^ (9.2.2) 

K, f 8 

g=l 

where S^^ Is Kronecker's delta; = 1 for g = 1, zero otherwise. The group 
diffusion equations are solved for cylindrical geometry with 0 and 3 symmetries 
by using a finite difference technique which is described in module RS-2. 
Fluxes in each group and each mesh-point are calculated by solving a system of 
linear equations. 



-c^(t,^ , + b^cf)^ - a^(t)^ = (9.2.3) 
n^n-1 n n n^n n ^ ' 



-^M*M-1 = ^M-l 



where indices 1, M denote mesh-points at the center and outer boundary of the 
reactor respectively. is the source at the point n. 

In matrix form the system (9.2.3) can be written as 
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(9.2.4) 



g O 

where ((> is a vector representing fluxes (j)^ at all mesh-points. The source 
term is a vector whose components are 



n n sn ^n K f n^n'' 



where h is volume of the mesh- interval. In matrix notation can be written 
n n 

as 



K, 



(9.2.5) 



where 



6. is the Kronecker delta and 



* = 



= (E*^^ E^^ 0) 

1 2 
F = (F-^ F F^) 



The vector ()>^ is formeo from the elements 4'^, the diagonal matrix E^^ is formed 

p— 1 a -1 

from the elements Z ' , and the diagonal matrix F*^ is formed from the elements 

sn 



j 



Finally, the assembly of group euqations may also be regarded as a matrix 
problem, which may be written as 



A 1 = s , 



(9.2.6) 



where 



A-" 0 0 




9,3 The Inner Iteration ,r^ ^ 

The inner iteration allows ca^'culation of spatial dependence of fluxes 
with respect to the boundary conditions. For simple one-dimensional geometries 
the /nner iteration can be avoided and fluxes at each mesh-point calculated by 
a direct procedure such as Thomas algorithm which is described in Module RS-2, 
Chapter 2. 'ihis algorithm is used in this code. 



9.4 The Outer Iteration 



The overall matrix problem to be solved is expressed in Eq. (9.2.6) as 



A i = S (9.4.1) 



in which 



S = E _$^+ i F.^^. - (9-4.2) 
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The normal procedure is to begin the solution of the group equations in 
the highest energy group by using a flm. guess for the calculation of the fission 
source. In the highest group there are no slowing-down scattering terms. Pro- 
ceeding to the next group, the slowing-down part of the source can be calculated 
from the new fluxes in group 1; since fission source in this model appears only 
in the highest group, the fission part of the source in group 2 is zero. The 
calculation can proceed in this fashion to thermal energy. At the end of the 
outer iteration we have a completely new flux vector, so that the fission source 
F may be recalculated and the whole process may be repeated. 

The equation which must be solved is 



(A - E)0 = ^ F ^. (9,^^ 3) 



The iteration process proceeds by defining a vector 



*=k^*=^^ (9.4.4) 



where 

S = F(},. 

(Note that S denotes the fission source only.) 
The iteration is defined by the relation 



(9.4.5) 



From Eq. (9.4.5) we can write 



* (9.4.7) 
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By combining Eqs. (9.4.5) and (1.4.7), the p-th iterate of the source 
reactor is given by 



(9.4.7) 



To estimate the eigenvalue one utilizes the fact that when a convergence is 
being approached, the fission source for one iteration should equal the fission 

source for the next. In this code the estimate of the eigenvalue is determined 

(0) (0) 
in the following fashion. The iteration starts with k equal to 1 and S 

such that 



S dv 



(vol of core) 



(9.4.8) 



Ther from (1.4.7) and (1.4.4), 



(1) 



S^^^dv. 



(vol of core) 



(9.4.9) 



and similarly 



,(p) _ 



(9.4.9) 



The iteration continues until convergence of the eigenvalue k^^^ is obtained, i. 



< e . 



(9.4.10) 
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The single power iteration described above tends to converge rather slowly 
and acceleration techniques are used to improve the convergence rate. If we 
denote the unaccelerated fission source for iteration p by S*^^\ then we 
calculate an accelerated fission source from 



S^P^ = S^P-^^ + u,(S*^P^ - S^P-^^. (9.4.11) 



This form of acceleration is referred to as a first order acceleration, Tt^e 
Chebyshev polynomial method is described Reactor Statics Module, RS-8, The 
optimum value of lo requires a knowledge of the maximum eigenvalue of the matrix 
(A - E)F Generally this eigenvalue is not known. One of the methods to 
solve this problem is described in the following section. 
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9.5 The Exponential Over-Relaxation Technique 

The program uses a special technique for the calculation of the source 
Introduced for the subsequent Iteration step. 

Fission sources are computed from the calculated fluxes at each mesh 
point as follows 



n J f^n n n = 1,2, . . . ^^.^.1; 



For the first Iteration the Initial source-guess will be, e.g.. 



(9.3.2, 



where R Is the radius of the core. The source for the p-th Iteration will be 
given by formula 

S^P^ = S^P-^^ + i4S* - S^P-^\ (9 5 3) 

n n n n ' 

where o) Is the relaxation factor whose optimal value lies between 1 and 2. 
The equation (9.5.3) can be rewritten as 



n 



S* 

n 



(9.5.4) 



n 
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The optimal value of the relaxation factor w is automatically computed 
during the calculation. 
Letting, 



n 



,(p-i) 



-1, 



(9.5.5) 



we can have two cases: 
1) Q' > 0. 

The exponential function can bo expressed as 



S 

exp[.(^!-^j -1] = 1 + ^(-^ -1) 
n n 



S* 

n 



(9.5.6) 



where 



Q = 03* 



,(p-l) 



-1 



(9.5.7) 



2) Q' < 0. 

In this case we have 



S* 

exp[a)(-^^— y -1)] = 



(9.5.8) 
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The calculation of the exponential function in the core is speeded up by the 
following approximations: 

For Q > 0.1 e^ exp[-Q] 



0.001 < Q < 0.1 -^ l + ln 



Q < 0.001 e^ 1 - Q 



In order to find the optimum value of o) it is necessary to scan the 
maximum relative error at each iteration given by the expression 

s(p) « s^P'^^ 

e^^\s) = max(-i-^p3j ) , (9.5.9) 

i 

where i runs through all mesh-points. The computation is normally started with 
03 equal to 1.9. At the end of 20 iterations the convergence ratio is tested. 
If the convergence is fast enough, let us say, if 

> . (9.5.10) 

the iteration with unchanged w is continued for another 20 iterations. If ^ 
the convergence is slow with respect to the above criterion, i.e., the inequality 
is reversed and if the relative point-wise error oscillates which can be 
described by the inequality 



1.0, (9.5.11) 
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w is reduced by 0.1 and with this new value of to the iteration process is 
continued in the same manner as described above. 
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9.6 Accuracy Criteria 

In this code two criteria for the outer Iteration are adopted. The 
convergence of the multiplication factor k will be reached when the following 
holds. 



^ 5.0 ' 



(9.6.1) 



(p-1) 



where k^*" is the eigenvalue obtained in the previous iteration. Once 
criticality is obtained, the program proceeds to obtain the pointwise con- 
vergence of neutron group fluxes which will be reached when the following 
condition is fulfilled. 



n n 



(p-1) 



(9.6.2) 



max 



With the poison search another convergence criterion must be Included which 
is discussed further. 
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9.7 Control Search 

The program does a control search on a poison or any fictitious material 
which adjusts the reactor to the critical condition. Let the multiplier of 
the poison concentration be denoted by X. In general, the parameter X can be 
applied to any other type of control. By changing X the program finds k, then 
changes X, finds k again, and so on. In the early stages of the process the k 
calculations stop when both criteria (9.6.1) and (9.6.2) are satisfied. Then 
the criterion to guarantee criticality is applied, i.e.. 



if k is chosen to be unity, the reactor is critical. Clearly k also can be 
c c 

chosen to take into account the presence of other reactivity controls • 

The initial guess for the multiplication factor proceeds according to the 
following formulas: 

a) If the fiist guess for the multiplication factor X is zero, the second 



I k(converged) - k | 



< G 



(9.7.1) 



guess Is 



X2 = AX *{of sign of (1.0 - k^)} . 



(9.7.2) 



b) If the first guess is non-zero the second guess is given by 



X2 = 1.0 + AX *{tsign of (1.0 - k^)]}* X, 



(9.7.2) 



where k is the converged eigenvalue and AX is an increment in the 
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parameter X, This formula tries to ensure that X moves in the right 
direction: it assumes that if X is negative, an increase of X 
decreases the reactivity, and the converse of this if X is positive 
(this would be applied, e.g., for the search of fissile isotopic 
composition). The programs finishes when three conditions (9.6,1), 
(9.6.2), and (9.7.1) are satisfied. 



With the quantities Xj^ and X2 the search is begun and X is varied in order 
to make k equal to unity. It is easy enough then to make a simple, search 
procedure such as ordinary "regula falsi," but these processes appear to be 
very inefficient in this context. The difficulty is that given X, we can 
find k but the more accurately we require it, the greater the number of itera- 
tions is required. One way to go is to improve "regula falsi" process so that 
the search is accelerated. In ordinary "regula falsi" at each stop, it is one 
of the limits (a, b) which is used and the last approximation of the root. 
Instead of that we can use the last two approximations which are Closer to the 
root which is searched for than the limits of the interval (a, b) . The formul 
which is used has ti..' form given by Equation (9.7.5) and is illustrated in 
Figure 1. 



a - a , 
n n-1 



n+1 n n f(a ) - f(a 



(9.7.4) 



n-1' 



is 



- 16 - 




The aj^ is found from the formula 



^1 = ^ - ^^^^ f(b? : f(a) • (9.7.5) 



Depending on the sign of f(a), f(b), and f(aj^), computed from the formula 

(9.7.6) if f(a) < 0, f(b) > 0 and fCa-^) < 0; if f(a^) > 0 from the formula (9.7.7), 



and 



b - 

^2 = ^1 - ^(^^ f(b) - f(aj (^•7-^> 



a^ - a 

^2 = ^1 - f(a,) - f(a) • (^-^-^^ 
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If perchance the point 33 will appear Out.ide the interval (a, b) then in the 
following step it will be necessary to consider the closer limit instead of 
the point a^. It can be shown that convergence of this method is faster than 
the ordinary "regular falsi". 
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9.8 Power Normalization of Neutron Fluxes 

The source-normalized power at the point r is given by the expression 

p(r) = ^^Kj4(?)*^(r) j = 1,2,3 (9.8.1) 

where j stands for energy produced per fission reaction, *^(r) is the source- 
normalized flux. In order to get fluxes normalized to the total power output 
of the reactor we have to multiply the source normalized fluxes by the normali- 
zation factor Y given by the formula 



Y = 




(9.8.2) 



0 

where dr = 27r rdr and H is the height of the reactor. The power-normalized 
flux is then given by 



J j 



(9.8.3) 
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9 . 9 Organization of the Program ^ 



The code WANDIC has the following limitation^; 
a) Number of energy groups is three 



b) The reactor can be divided into four regions from which the outer 
one is non-multiplicative. If use of fewer regions is desired, the 
the entry data for some of these regions will be identical. 

c) Memory array space is provided for lip f9j30Qvmesh--points. 

d) Number of intervals dividing each region must be even. 

The code consists of a number of subprograms which perform different parts of 
calculation and are controlled by the main program. 



9.10 The Role of Different Subprograms ^ ^ \: 

PLOT: The subprogram participating ip preparation of the output and 
plotting the related results on the printer. 

SIMP: The subprogram performing the integration 'of functions over the 

^> 

core volum. by SIMQSON'S rule. 

ZERO: The subprogram which employs the accelerate%i;fi:egula falsi" tech- 
nique. •J*' 



9.11 Flow-Chart ' " - 

Besides its role as the manager of subprograms, the main program is in 
charge of entry of data for calculation of fluxes and preparation of results 
for the output. In Figure 2 the f low-chat|r ,of ^the sequence of calculation is 
presented. ef/ 
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FLOW-CHART A 




INITIAL SOURCE 
S(I.N). S(2.N). S(3.N) 



kant = l 



DO 
I = 1.2.3 



CALCULATION OF 
<I>(I,N) 



< 



CALCULATION OF 
S(2,N) 



> 





2 


CALCULATION OF 
S(3,N) 







■c 



CALCULATION OF 
S(I,N) 



DO END 



YES NO 



> 



k=/s(r) dr 



"EXPONENTIAL 
OVERRELAXATION" 



ay 



NO 



EIGENVALUE 
CONVERGENCE 



YES 



ITER- 



NO 



FLUX 
CONVERGENCE 




FIG, 2 CON'T. 
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6ENERAL FLOW-CHART OF THE PROGRAM 



NIT=4 



ZERO 
("REGULA FALSI"; 



NIT : 1.2, 



3E 



\ NO 



(^Ak»(Ak)Qn^-o)-^ 



Ak-(Ak)an^- 



NIT=3 



FIG. 2 




INPUT 
READING OF 
DATA 



ITER =1 
NIT = I 
X = I 

(AN INITIAL VALUE 
FOR A CHOSEN 
PARAMETER IS SET) 



NIT=2 



DISTRIBUTION OF 
MESHPOINTS. 
CALCULATION OF 
MESH-SIZE H(K) 



X«X-AX 
1 



CALCULATION OF 
FLUXES. INNER ITE- 
RATION. CALCULATION 
OF k. (FLOW-CHART A) 



=^RINT k, X 



^ Ak \ 

TEST FOR CRIT.CALITY 
V CONVERGENCE / 



c 



ITER: 







CALCULATION OF 
POWER NORMALIZED 
FLUXES 






CALCULATION OF 
OUTPUT RESULTS 






PLOT 
(PRINTING OF 
RESULTS) 






OTHER CASES? ^ 



> 



Ak: 



3 



ITER "2 



YES 



JNO 
STOP ) 
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12 Output 



The output of the program are: 

a) source normalized group fluxes versus radius 

b) relative magnitudes of fluxes at each mesh-point with respect to the 
flux at the S3rmmetry axis of the core 

c) ratio of group fluxes 

d) eigenvalue of the unpoisoned reactor 

e) poison concentration needed for criticality 

f) power normalized group fluxes 

g) power generated in each fueled region with respective power densities 

h) for user who wants to watch the iteration process, printing of eigen- 
value K computed at each iteration is provided. 

13 i.ames of Principal Variables Used by the Program 
mesh-size of the region K 

neuLron flux of the group I at the mesh-point N 
slowin^-down source at the mesh-point N 
fission source at the mesh-point N 
source at the mesh-point N 

maximun relative error for the source on the p-th iteration 

K = 1 
c 

multiplication factor for the control poison cross-section 
eigenvalue of the preceding iteration 
power density of the region K. 



H(K): 

FLUX(I,N): 

SDS (N) : 

FS(N): 

S(N): 

EPS(P): 

FME: 

XP: 

FMANT: 
FPR(K) : 
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^•l^ Use of the Program for a Simple ProbJam, Presentation of Entry of Da ta 

The data Introduced here are for Illustration purposes only. Consider 
a cylindrical reactor consisting of three different fissionable regions and a 
reflector (see figure below). The outer radius of the first region Is = 87,7 
the second R^ = 124,00 cm, R, = 152,00 cm, R. = 159 cm. 



R = 0 


REGION 
I 


REGION 
II 


REGION 
III 


REGION 
IV 




R^ = 0 


^2 


h 




^5 





The data for the problem would start: 
1 card: symbols used for reading are 

OVREL CONV PREC BH2 SCP POW NINT(l) NINT(2) NINT(3) 

NINT(4) ICP ICC 

where 

OVREL: the relaxation factor (Sec, 1*5). This factor for a new problem 
is put equal to 1.9. Once an "optimum" factor is found, this 
may be introduced into a set of similar problems (e.g,, with the 
same mesh-point distribution), 

CONV: the pointwise flux convergence accuracy number (e«), (See Eq, 
(9*6,2) ^ 

PREC: the accuracy (c^) for the eigenvalue convergence, (See Eq, (9,6,1) 
BH2: axial buckling, 

SCP: microscopic cross-section of a poison material for the third group 
energy level in barns. 

POW: thermal output of the reactor in watts. 
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NINT(K): number of intervals per region K. 

ICP: control search parameter. If its value is 

1... control poison is considered in region 1 

2. • .region 2 

3. . . region 3 

4... all fueled regions 

5... regions 1 and 2 

6... regions 2 and 3 
ICC: is the case control parameter (equal to zero, if only one case 

is to be treated; positive, if otherwise) 

The three group diffusion equations in one of the multiplicative powers 
assumes, in general, the following form 



2 3 (vZj,*, 



V'*3 - ^a3 *3 ^s,2.3*2 = ^ 



For example, with calculatad macroscopic cross-sections 

1.7 A^*^ - (0.004 + 0.06)*^ + i (0.43 *^ + 0.015 + .16 *^ = 0 
0.85 A^$2 - (0.025 + 0.05)*2 + 0.06 *^ = 0 
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0.39 th^ - 0.13 *2 + 0.05 *2 



Let us divide the region I into 30 equal intervals, the regions II and III into 
20 intervals and region IV into 4 intervals. The mesh-points in the radial 
direction will be numbered as showti in the figure below. 

1 31 51 71 75 

4 ^ \ 1- A \ ^ 



R = 0 « 87.7 cm R^ = 124.00 cm R, = 152.0 cm Re = 159.0 cm 

1 2 3 4 5 



-3 2 

Let us assume that the axial buckling BH2 = 1.56 * 10 cm ; we shall use poison 
search in all fissionable regions, ICP = 4, the microscopic absorption cross- 
section of the poison material a = 2109 bn, and this will be the only case to 
be treated. The power output of the reactor is 2441 MW. 

OVREL CONV PREC BH2 SCP POW 

1.9E - 01 l.E - 02 l.E - 04 156. E - 05 2109. E 00 2441. E + 06 



ICP ICC 

30 20 20 04 04 0 



2, 3, 4th Cards : 



The second card carries physical data for the first region — group 1, the 

third those for group 2, the fourth for group 3. Data are read in the following 

wa t ts 

order (the symbols are self-explanatory; means fj^ggj^Q^^ 
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2nd: 17. E - 1 4.E - 3 6.E - 2 43. E - 2 93. E - 15 

3rd: 85. E - 2 25. E - 3 5.E - 2 l.E - 2 19. E - 14 

4th: 39. E - 2 13. E - 2 O.E + 0 16. E - 2 21. E - 13 

Other 6 cards carry the similar data for regions II and III. 

11, 12. 13th Cards ; 

These cards carry the data for the outer non-fissionable region in order 
of groups 1, 2, 3. These data are 



^aj h j-^j+l 



Let the reflector have the following data in order of groups 1, 2, 3 

11th: 19. E - 1 36. E - 5 26. E - 3 
12th: 92. E - 2 57. E - 5 28. E - 3 
13th: 31.^ - 2 85. E - 4 O.E 0 

14th Card ; 

This card carries outer radii of regions with increasing order 



*2 h h »5 

877. E - 1 124. E 0 152. E 0 159. E 0 
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15th Card ; 

The value carried by this card is the increment of the multiplier adjusting 
the poison concentration at early stages of control search. 

PP 
0.05 



9.15 Rules for Card Punching 

The format of data read is fixed: real variables are read with the format 
E 10.0, the integer ones with the format 12 with one blank space separating them. 
On the last card PP^ variable is read with the format F5,3. 
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MAIN DATE = 74327 



12/56/19 



0001 

0002 
0003 
0004 
0005 



) 



0006 
0007 



0008 
0009 



I 
C 
C 
C 



HA IN PROGRAM 
PROGRAM HANDIC 



THE PROGRAM SOLVES THE THREE-GROUP DIFFUSION EOUATICIfo 

FOR A CyLlNORICAL REACTOR CONSISTING CF THREE FUELED RcGlUNS 

AND A REFLECTOR 



INPUT-MEANING Of SYMBOLS 

NINT... NUMBER OF INTERVALS FOR EACH REGI.IN 
OVREL...OVERRELAXATION FACTOR 

CONV. ..PRECISION FOR THE POINTWISE FLUX CONVERGENCE 

X... MULTIPLIER FOR THE CONTROL SEARCH PARAMETER 

0.. .DIFFUSION COEFFICIENT 

SCA., .MACROSCOPIC ABSORPTION CROSS-SECTION 

XNF... MACROSCOPIC FISSION CROSS-SECTION+NU 

XKF.. .MACROSCOPIC FISSION CROSS-SECTION*KAPPA 

SR. ..SCATTERING CPQSS-SECTIQN INTO THE LOWER ADJACENT ^KCIUP 

BH2.. .TRANSVERSE SUCKLING 

RR.. .OUTER RADIUS OF THE REGIONS l,2,3f't 

POrt.. .POWER OUTPUT OF THE REACTOR 

SCP.. .MICROSCOPIC CROSS-SECTION OF THE CGNTr^OL POISJN 

ICP...C0TROL PARAMETER FOR THE POISON SEARCH 

1... REGION 1 

2. ..REGION 2 

3. ..REGION 3 

i»,..ALL FUELED REGIONS 

5. ..REGIONS 1E2 

6. ..REGIONS 2C3 

ICC. ..PROGRAM CONTROL INDICATING THE NUMBER OF CASES 
POSITIVE INTEGER. ..SEVERAL CASES TO BE TREATED 
O...ONLV ONE CASE 



LIMITATIONS OF THE CODE 

1ST:NINT(KI MUST BE EVEN NUMBER 

2ND: TOTAL NUMBER OF MESH-POINTS BE LESS THAN 301 

SUBROUTINE PLOT (LI.N.LP) 

REAL MX,HNiPRNT{6,9),DISP(6l,DY,X(300)fY(3,300),H(4) 
INTEGER LI,LP,PT,STRG(10n,BLK,PER(21).AST,I,N,NINTIA) 
COMMON y.HiNINT 

DATA PRNT/4HFAST,4H FLU,4HX - ,4HGR0U,4HP 1 ♦ IH ,4HRt^'jU,4HNANC,4HE 
1 FL»4HUX -,4HGR0U,4HP 2 ,4HSL0H,4H FLU,4HX - ,4HGR'' ',4HP 3 ,IH ,4H 
26R0U,4HP 1 ,1H ,1H ,1H t IH ,4HGR0U,4HP 2 ,1H ,1H ,LH ,1H ,4HGR0U.4 
3HP 3 flH ,1H ,1H tlH ,4HGRCU,4HP I ,4H/ GR,4H0UP ,4H3 t IH ,4HGR0 
4Ut4HP 2 .4H/ GRi4H0UP t4H3 »1H ,4HGR0U,4HP 1 ,4H/ GR,4H0UP ,4H2 
5 ilH /iaLK,AST/lH , IH*/, PER/21*' . •/ 

DO 1 1=1,101 

STRGin=6LK 

CALCULATION OF RADII CORRESPONDING TO THE MESH-POINTS 



A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 

A 
A 
A 
A 
A 

:^ 

A 
A 
A 



XID^O.O 
K=l 



3\ 



A 
A 
A 
& 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 



I 
2 
3 
4 
5 
6 
7 

a 

9 

10 

11 

12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
,27 
28 
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30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 
51 
52 
53 
54 
55 
56 
57 
58 
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I 



0010 




DO 2 I=2tN 




A 


59 


0011 




X(I)=X(I-11*H(K) 




A 


60 


0012 




IF (I.EO.NINT(K)) 


K=K + 1 


A 


61 


0013 


2 


CONTINUE 




A 


62 


0014 




WRITE (6t9) (PRNT{I,LP)fI=li6) 


A 


63 


0015 




MN=Y<LI,1) 




A 


64 


0016 




MX-MH 




A 


65 


0017 




DO 6 I=2tN 




A 


66 


0018 




IF (MX-YlLIfD) 3 


,4,4 


A 


67 


0019 


3 


MX=Y(LI,n 




A 


68 


0020 




GO TO 6 




A 


69 


0021 


4 


IF (MN-Y(Llfl)) 6 


,6,5 


A 


70 


0022 


5 


MN=Y(LI,I) 




A 


71 


0023 


6 


CONTINUE 




A 


72 


002^ 




OYsiCX-MNi/lOO 




A 


73 


0025 




DO 7 1=1,6 




A 


74 


0026 


7 


DISPII)=MN+(I-1)*20*DY 


A 


75 


0027 




WRITE (6flO) DISP 




A 


76 


0028 




WRITE 16, U) (PER(niI=li21J 


A 


77 


0029 


; 


00 8 1 = 1, N 




A 


78 


0030 




PT=(Y(LI,I)-HN1/0Y+1.5 


A 


79 


0031 




STRG(PT)=AST 




A 


80 


0032 




WRITE 16,12} XII) 


,5TRG 


A 


81 


0033 


8 


STRG{PT)=BLK 




A 


82 


0334 




RETURN 




A 


83 




c 






A 


84 


0J35 


9 


FORMAT {1H1,45X,7HPL0T 0F,IX,6A4,//,57X,6HZ-aXISI 


A 


85 


0036 


10 


FORMAT (11X,6(E10. 3,10X1) 


A 


J*6 


0OJ7 


11 


FORMAT I16X,2IIA1 


,4X1 ,/,7H R-4XIS1 


A 


87 


C033 


12 


FORMAT (fcll.3,5X, 


lOlAl) 


A 


88 


0039 




END 




A 


89 



32 
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c 


INTEGRATION SIMPSON FOR N POINTS. N MUST BE EVEN* 


Q 

D 


1 
i 




c 


N MUST WITHIN 4 AND 100 


n 


-> 


OOOL 




FUNCTION SIMP (DXiXXiNtJl) 


Q 




0002 




OIHtNSlON XX(100)t DX(4) 


n 
D 




0003 




SIPPj^XXU )+XX(N)*4«*XX(N-ll 


D 
D 




000* 




J={N-3)/2 


Q 

U 


0 


0005 




00 1 I^ltJ 


B 


7 


0006 


1 


SIMP=SIMP+^.*XX(2*I )+2.»XX(2*I*l 1 


n 

P 


A 

0 


0007 




SIMP=SIMP*0X(Jl)/3. 


B 


9 


OOOB 




RETURN 


B 


10 


0009 




END 


B 


11 



33 



ERIC 
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C 


ACCELERATED 'REGULA FALSI' 


C 


1 




c 




C 


2 




c 




c 


3 


0001 




SUBflOUyiNE ZERO (Al,A2,FlfF2f AlO,A20,DFl,DF2rRAI)) 


c 


4 


0002 


I 


RAD=A2-F2*(A2-Al)/(F2-Fl) 


c 


5 


0003 




VARA*A2-F2*(A2-Al)/IF2-Fi) 


c 


6 


000^ 




HH=ABSCA20-fllO) 


c 


7 


0005 




Hl=ABS( VAKA-AIO) 


c 


8 


0006 




H2=ABS(VARA-A20) 


c 


9 


0007 




IF IH1-H2) 2f8,3 


c 


10 


0008 


2 


TEMP=H2 


r 


11 


0009 




L=2 


c 


12 


0010 




60 TO k 


c 


13 


0011 


3 


TEMP=Hl 


c 


14 


0012 




L=l 


c 


15 


0013 


<t 


IF ITEMP-HH) 8i5,5 


c 


16 


OOU 


5 


GO TO (6,7)t L 


c 


17 


0015 


6 


A2=A20 


c 


18 


0016 




F2=OF2 


c 


19 


0017 




GO TO 1 


c 


20 


0018 


7 


A2=A10 


c 


21 


0019 




F2=DF1 


c 


22 


0020 




GO TO 1 


c 


23 


0021 


8 


RETURN 


c 


24 


0022 




END 


c 


25 



34 



ERIC 
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' C ENTRY QF DATA AND DIMENSIONS OF ARRAYS D I 

0001 DIMENSION NINT(<t)f RRI5], H(4)f EPSI30), SCAOI3,3r D 2 

I 0002 DIMENSION SCA(3,4), SRI3,4), D(3fA), XNF(3,3), XKF(3,3I 0 3 

' 0003 DIMENSION SANTI300), DFMAX(3), FSI300), SUM(300)i SDS(300), AIDD(3 D 4 

If 300) D 5 

% 0004 DIMENSION ALFA(300), BETAOOO), S(300)i FLUX(3,300), (iXdOO), FPR( 0 6 

' 13) D 7 

0005 COMMON AIDD,H,NINT 0 8 

% OOOb 1 READ (5.142) OVREL ,C0NV,PREC,BH2, SCP.POWi NINT , ICP ,ICC D 9 

' 0007 READ 15,143) ( ( 0( Ii K) jSCAU i K ) ,SR{ h K), XNF( I,K) ,XKF( I ,K) , 1 = 1 ,3) ,K= D 10 

11,3) 0 II 

I 0008 READ (5,144) ( D( 1 .4) , SCA( 1 ,4) ,SR( 1,4) tl = l,3) D 12 

0009 READ (5,145J (RR(K) ,K=2,5) 0 13 

0013 READ (5,I46J PP 0 14 

I *C D 15 

C PRINTING CF DATA D 16 

C D 17 

\ 0011 RR(l)=0,0 0 18 

0012 WRITE (6,147) D 19 

0013 WRITE (6,148) >J 20 
) 0014 WRITE (6,149) (RR(L ) t L= 1,5) D 21 

0015 WRITE (6,150) NINT,CVREL,PREC,CCNV D 22 

0016 WRITE (6,151) BH2,SCP D 23 
) 0017 WRITE (6,152) ( (SCA( I ,K) , I=1,3),K=1,4) D 24 

0018 WRITE (6,153) ( (D( I ,K ) ,1=1,3) ,K=l ,4) 0 25 

0019 WRITE (6,154) ( (SR( I ,K) , 1 = 1 ,3 ) ,K=l,4 ) D 26 
) 0J20 WHITE (6,155) ( (XNF( I ,K ) , 1=1 ,3) ,K=1,3 ) 0 27 

C02I WRITE (6,156) ! ( XKF( I ,K) , 1=1 ,3) ,K=1 ,3 ) D 28 

0J22 kRITE (6,157) POW D 29 

I 0023 WRITE (6,158) H 30 

C 0 31 

C NUMBERING OF MESH-POINTS D 32 

) C 0 33 

0024 PI = 3. 14159265 D 34 

0025 CF=1.6021E-13 0 35 
) 0026 TEMP=1.0 0 36 

0027 ITER=1 D 37 

0028 NIT=1 0 38 
) 0029 MNT(1)=N1NT(1)+1 0 39 

0030 N1=NINT(1) D *0 

0031 00 2 K=2,4 0 41 
) 0032 2 NINT(K)=NINT(K)+NINT(K-1} 0 42 

0033 N2=NINT(2) 0 *3 

0034 N3=NINT(3) 0 

% 0035 N3M1=N3-1 0 *5 

0036 N4=NINT(4) or ^ 

0037 N4M1=N4-1 OD ^7 
I 0038 XP=1. 0 ^8 

0039 JJ=1 0 

0040 DO 3 K=l,4 0 50 
J 0041 00 3 1=1,3 0 51 
' 0042 3 SCA(I,K)=SCA(I,K)fSR(l,K)+BH2*D(I,K) D 52 

0043 DO 4 K=l,3 D 53 

X 0044 4 SCAa(3,K)=SCA(3,K) D 54 

' 0045 DO 5 N=l,N4Ml D 55 

0046 FS(N)=0.0 0 56 

^ 0047 5 SUM(N)=0.0 0 57 

ERIC t D 58 
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0048 
0049 
0030 
0051 
0052 
0053 



0054 
0055 
0056 
0057 
0058 
0059 
0060 
0061 
0062 
0063 
0064 
C065 
0066 
0067 



0J6d 
0069 
0073 
0071 
0072 
0073 
0074 
0075 
0076 
0077 
C078 
0079 
0080 



ERIC 



0081 
0082 
0083 
0084 
0085 
0086 
0087 
0088 
0089 
0090 
0091 
00<)2 
0093 



7 
8 
C 
C 
C 
9 
10 

u 



12 



13 

14 
15 

C 

C 

c 

16 



17 
Id 

19 
20 
21 

C 
C 

C 

C 

22 



23 



DO 8 K=l,4 

IF (K-1) 8,6,7 . • 
TEHP=NINT(ll-l ■ /' 

GO TO 3 ' 

TEHP=NlNTJ<<r-MNT(K-l) /' 
H(K) = 1RR(K*U-RR(KJ)/TEKP 

POISON CONTROL BLOCK 

GO TO (10,10,10,11,12,13}, ICP 

SCA<3,ICP)=SCA0(3,ICP)»XP 

GO TO 16 

KF=l 

KL=3 

GO TO 14 

KF=1 

KL=2 

GO TO 14 
KF=2 

KL=3 c-r 

DO 15 ^=KF,KL 
SCA(3,K) = SCA0(3rK)>»XP 
JJ = 1 



A. 



GE^Ji:RATION OF THE INITIALtFAST JOUPCF' 

S(ll = l./(R»^(4)*PK('»)'PL)Jy: ,, 
ITC=1 ' ■'■'v 

SANT(1)=S(1) 
DO 21 N=2,N4 
IF (N-N3) 18,19,17 
SIN)=0.0 
GO TO 21 
S(.\I = S(1J 
GO TO 20 

S(N3)=S(l)*H(3)/(h(3l+H(4l) 
SANT<N)=SIN> 
CONTINUE 
FH4NT=1.0 



CALCULATION OF FLUXES 

00 58 1=1,3 r 
TEMP=H(l)»H(l)*SCA(I,l)H.*D{I,l)i: 
ALFA(2)=4.»D(I,1)/TEMP ' 
B ET A 1 2 ) =S { 1 1 *H m *H m / T EfJP ^-^ 
R=0.0 g 




36 



DO 26 N=2,N4M1 
IF (N-NINT(K)) 23,24,23 
R=R4H(K) 

AN=(l.+H(K)/(2.*R)J*0<l,K)/(H(K,)*H(K U 
CN=(1.-H(K)/(2.»R))*0<1,K)/(H<K)«H(*>1 
BN=AN*CN+SCA(I,K) , 

GO TO 25 , .* I 



D 
0 
0 
D 
0 
D 



59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 
70 
71 
72 
73 
74 
75 
76 
77 
78 
79 
80 
81 
82 
83 
84 
85 
86 
87 
88 
89 
90 
91 
92 
03 
94 
95 
96 
97 
98 
99 
100 
101 
102 
D 103 
0 104 
D 105 
D 106 
0 107 
D 108 
0 109 

n no 
111 

112 
113 
114 
115 
116 
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0094 
0095 
0096 
0097 
0098 
0099 
0100 
0101 
0102 
0103 
0104 
0105 
0106 
0107 
0108 
0109 
Oil) 
0111 
0112 
0113 
0114 
0U5 
0116 
0117 



24 K=K*1 
R=RR(K) 

TEHPs(H(K-l)+H(K))*.5 
AN=(1.+H(KJ/(2.»R))»D(I,K)/(H|K)*TEHP) 
CN=(l.-H(K-l)/(2.*RJI*0(I,K-ll/(H(K-ll*TEMPI 
BN=ANtCN*(SCA(IiKI*H(K)*SCA(I,K-n*H(K-l))/«H(K|*H(K-lJ) 

25 TEMP=BN-CN«ALFA(N) 
ALFA(N+1)=AN/TEHP 

26 BETA{N+1)=(S(N)*CN»BETA(N)»/TEMP 
BM=CN»l3.-»H<4)/(l.066*D(Ii4))l-AN 
CM=4.*CN-BN 

GO TO (28,271, ITER 

27 TEMP=FLUX(I,N4) 

28 FLUX(I,N4)=(SIN4M1)*CH*BETA(N4))/(BH-CH*ALFA(N4)» 
GO TO (30,29), ITER 

29 0FHAX(n=ABS((TEMP-FLUX(I,N4>)/TEKP) 

30 00 35 N=1,N4M1 
L»N4-N 

GO TO (32,31), ITER 

31 TEHP=FLUX(I,L) 

32 FLUXII,L) = ALFA(L+l)*FLUX(I,'.+I) + eETA(L + l) 
GO TO (35,33), ITFR 

33 OFMAX(n=ABS((TEMP-FLUX(I,L))/TEMP) 
IF (DFHAX(I)-TEHP) 34,35,35 





fi 1 1 


34 


DFMAXd )=TEMP 






35 
r 


CQNTINJE 






r 
r 


CALCULATION CF SOURCE FOR 1 


) 


0120 




K=:l 




0121 




DO 58 N=liN4nl 




0122 




IF (N.GT.MNT(3)) GO TO 40 




0123 




IF (N-NINTIK)) 36,40,36 




0124 


36 


GO TC (38,38,37), I 




0125 


37 


SOS(N)=0.0 




0126 




GO TC 39 




0127 


38 


SOS(K)=SR(l,K)*FLUX(I,N) 




0128 


39 


FS(N)=XNF(I,K)*FLUX(I,N) 




0129 




GO TO 53 




0130 


40 


K=Kfl 




0131 




IF (K-4) 45,41,50 




0132 


41 


GO TO (43,43,42), I 




0133 


42 


S0S(N)=0.C 




0134 




TEMP1=0.0 




C135 




GO TO 44 




0136 


43 


SDS(N)=SR(I,K)*FLUX(1,N) 




0137 




TEMP1=SR(I,K-1)*FLUX(I,N) 




0138 


44 


TEMP2=XNF(I,K-l)*FLUXlI,N) 




0139 




FS(N)=0.0 




0140 




GO TO 49 




0141 


45 


GO TO (47,47,46), I 




0142 


46 


SOS(N)=0.0 




0143 




TEMP1=0.0 




0144 




GO TO 48 




0145 


47 


SDS{N)=SR(I,K)*FLUX(l,N) 




0146 




TEHP1=SR(I,K-1)«FLUX(I,N) 




0147 


48 


FS(N)=XNF(I,K)*FLUX(I,N) 


ERIC 


0146 




T£MP2'XNF(I,K-1)*FLUX(I,N) 



3 



i 



D 
0 
0 
D 
D 



0 
0 
0 
0 



0 
D 



n 
D 
D 
0 



0 117 
D 118 
D 119 
C 120 
D 121 
122 
123 
124 
125 
126 
0 127 
0 128 
129 
130 
131 
132 
0 133 
0 134 
D 135 
136 
137 
D 138 
D 139 
140 
141 
142 
143 
D 144 
D 145 
D 146 
0 147 
0 148 
U 149 
D 150^ 
D 151 
0 152 
D 153 
0 154 
D 155 
D 156 
D 157 
0 158 
D 159 
D 160 
D 161 
0 162 
163 
164 
165 
166 
167 
0 168 
0 169 
D 170 
171 
172 
173 



0 
D 
0 
0 
D 



0 
0 
D 



0 174 
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0U9 ^9 TEMP=H(K)+H(K-l) ^ 1^5 

0150 SOS<N)=(S0S(NI»H(KI+TEHPl*H(K-l))/TEHP D 176 

0151 FS(N) = {FS(N)*H(K) + TEMP2*H(K-U)/TEMP D 177 
C152 GO TO 53 0 178 
0153 50 GG TO {52i52.51), I P H"? 
015^ 51 SDS(NI=0.0 0 180 

0155 GO TO 54 D 181 

0156 52 SDS(NJ=SR(I,4)*fLUX(IiN) D 182 

0157 GO TC 54 f 183 

0158 53 SUH(N)=SJM(N)+FS(N) R 18^ 

0159 54 GO TC (56t56,55)i I D 185 

0160 55 IF (N.GT.NINT(3)) GO TO 57 D 186 

0161 S(NI=SJHIN) ^ 187 

0162 SUM(N)=0.0 D 188 

0163 60 TC 58 0 189 

0164 56 S{N)=S0SIN) D 190 

0165 GO TO 58 0 191 
0165 57 S(N)=0.0 D 192 

0167 58 CONTINUE ^ 1^3 
C 0 194 
C INTEGkATION of THt SOURCE D 195 
C D 196 

0168 FMULT=S(N3)»RR(4)*H(4)/3. D 197 
C169 00 63 K=l,3 D 198 

0170 GO TO (59,60t60), k 0 199 

0171 59 L=0 0 200 

0172 GO TC 61 0 201 
60 L=NINT(K-1)-1 0 202 

0174 61 M=MNT<K)-L ^ 203 

0175 R=RR(K) D 204 

0176 00 62 N=.l,M ^ 205 

0177 ITEMP=NtL 0 206 

0178 GX<N)=S(ITE^'PI*R 207 

0179 62 R=R+H(K) D 208 
C180 63 FMULT=FHJLT+SIMP(h,GX,H,K) 0 209 

0181 FHULT=FHULT*2.*Pl ^ 210 

0182 WRITE (6iI59) FHULT 0 211 

0183 IF (ABS(FfULT-FMANT)/FHULT-PRLC/5.) b5t64,64 D 212 
C 0 213 
C PREPARATICN OF ITEKATIONS 0 214 



C 



215 



0184 64 RPSP=0.0 ^ 216 

0185 DO 66 N=1,N3 0 217 

0186 S(N)=S(N)/FMULT D 2'8 

0187 EHAX=A(JS((S(N)-SANT(Nn/SANT(Nn 0^ t\ ^ 

0188 IF (EMAX-RPSP) 66,66,65 O ! /X- '^220 

0189 65 RPSP=EMAX ' ^ D 221 

0190 66 CONTINUE ° 222 

0191 EPS(JJ)=RPSP D 223 

0192 IF UJ-20) 72,67,67 0 224 

0193 67 IF UBS(EPSil0)/EPS{20))-l. 5) 68,71,71 0 225 

0194 68 IF (ABS(EPS(20)/EPS(19)*(l.UBS(EPS(19)n>- 71,69,69 D 226 

0195 69 GO TO (70,71), ITER 0 227 

0196 70 OVREL=OVRFL-0.1 D 228 

0197 71 JJ=1 0 229 

0198 72 DC 84 Nsl,N3 D 230 
cnYr 0199 Q»S(N)/SANHN)-1. D 231 
^ly^ 0200 IF (Q) 73,74,74 ^ 232 

\ 
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0201 


73 


Kl=l 


0202 




GO TO 75 


0203 


74 


Kl=2 


0204 


75 


0=A6S(0)*CVREL 


0205 




IF (0-0.1) 77,76,76 


0206 


76 


Q=1.0/EXPiO} 


0207 




GO TO 80 


0208 


77 


IF IQ-.OOl) 79,78,78 


0209 


78 


0=11. -0/2. l/n.+Q/2.) 


0210 




GO TO 80 


0211 


79 


Q=l.-Q 


0212 


dO 


GO TO (81,82), Kl 


0213 


81 


S{N)=SANT(N)*Q 


02 H 




GC TO 33 


0215 


82 


S(N»=SANT(N)*(2.-Q) 


0216 


83 


SANT(NJ=S(N) 


0217 


84 


CONTINUE 


0218 




FMANT=FHULT 


0219 




JJ=JJ+1 


0220 




ITC=ITC*1 


0221 




IF {ITC.CT.400) GO TO 140 


C222 




GO TO 22 


022J 


65 


GO TC (39,36), ITER 


0224 


86 


DFMM=0.0 


0225 




DO 68 1=1,3 


0226 




IF (DFMAX(n-OfMM) 68,86,87 


0227 


d7 


nFMM=OFMAX(I) 


022." 


dS 


CCNTINJE 


0229 




IF (DF"(M-CONV) 89,64,64 


0233 


89 


FME=1. 


0231 




DFMUL=FHULT-FHE 


0232 




V.RITE (6,160) XP.FMULT 


0233 




IF (XP.EO.1.0) E1GEN=FMULT 


0234 




IF {ABS(OFHUL)-PREC) 104,90,90 


0235 


90 


GO TC (91,92,99,103), NIT 


0236 


91 


0FM1=DFMUL 


0237 




NIT=2 


0238 




GO TO 93 


0239 


92 


IF (0FMUL*DFM1) 95,104,93 


0243 


93 


IF (OFMUL) 94,104,94 


0241 


94 


X1=XP 


0242 




XP=XP+PP*CFMUL/AbS(OFMUL) 


0243 




0FM1=DFMUL 


0244 




60 TC 9 


0245 


95 


IF (DFHUL) 96,104,97 


0246 


96 


DF2=DFM1 


0247 




DF1=0FMUL 


0248 




A20=X1 


0249 




A10=XP 


0250 




XP=X1-0FHI«(XI-XP)/(DFM1-DFMUL) 


0251 




GO TO 98 


0252 


97 


DF2=DFHUL 


0253 




DF1=0FM1 


0254 




A20=XP 


0253 




A10=X1 


0256 




XP=XP-DFMUL*(XP-X1)/(DFHUL-0FKI) 


0257 


98 


NIT=3 


0258 




K1=:XP 



38 



0 233 
0 234 
D 235 
0 236 
D 237 
D 238 
0 239 
0 240 
0 241 
0 242 
D 243 
0 244 
0 245 
D 246 
0 247 
0 248 
D 249 
D 250 
D 251 
[) 252 
0 253 
0 254 
0 255 
0 256 
D 257 
0 258 
0 259 
0 260 
D 261 
0 262 
D 263 
0 264 
0 265 
0 266 
D 267 
0 268 
0 269 
D 270 
n 271 
D 272 
0 273 



D 275 
D 276 
D 277 
0 278 
D 279 
P 200 
0 281 
0 282 
D 283 
n 284 
0 285 
0 286 
D 287 



D 289 
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0259 




ofmi=dfmul 


D 290 


0260 




GO TC 9 


D 291 


0261 


99 


NIT«4 


D 292 


0262 




dfhi=dfhul 


0 293 


0263 




Xl=XP 


0 294 


0264 




IF (DFMUL) 10Cil04,101 


D 295 


0265 


100 


XP=Xi-DFMl*(A20-Xl)/(0F2-DFMl) 


0 296 


0266 




GO in i02 


0 297 


0267 


101 


Xf-Xl-DFMlXXl-AlO/IDFMl-DFll 


D 298 


0268 


102 


XANT=XP 


0 299 


0269 




GO TO 9 


■ 0 300 


0270 


103 


CALL ZeRO {XUXANT,DFMlf0FMUL,Al0,A20»DFl,DF2,XP) 


0 301 


0271 




X1=XANT 


D 302 


0272 




XANT=Xi 


0 303 


0273 




0FM1=DFMUL 


D 304 


0274 




GO TC 9 m 


D 305 


0275 


104 


GO TO (105,106) , ITER 


D 306 


0276 


105 


ITER=2 


0 307 


0277 




WRITfc |6,16U ITER 


D 308 


0278 


C 


GO TO 6^ 


D 309 
0 310 




C 


PREPARATION CF OUTPUT 


0 311 




C 


CALCULATION OF POtJEK NORMALIZED FLUXES 


0 312 




C 




D 313 


0279 


106 


FP=0.0 


D 314 


0280 




DO 107 Nsl,N3 


n 315 


0^31 


107 


SUM(NjrfO,0 


n 316 


028J 




K=l 


D 317 


0283 




DO 110 I=l»3 


D 318 


0234 




DO no N^liNSf*! 

IF INS<lNT{Kn lC8fl09,l08 


D 319 


0285 




D 320 


0286 


106 


FS(N)=XKF(IfK)*FLUX{I,N) 


D 321 


0287 




GO TO HO 


D 322 


0288 


109 


K=K+1 


0 323 


0289 




TEfP=)tKP( I,K-ll*FLUXl I,N) 


D 324 


0290 




FS(N)=)IKF<I.K)*FLUX(I,N) 


0 325 


0291 




FS(N)*|f=S(N)*H(Kl+TEMP»H(K-l))/l"'K)tH(K-l)) 


0 lib 


0292 


110 


SUMINMSUM(N)*FS(NI 


0 327 


0293 




TEMP=h(3J/(H(3)*H(4)) 


0 328 


0294 




DO m I=lf3 


0 329 


0295 


111 


SUM(N3)^SUMJN3J*XKF( I ,3)*FLUX( 1 ,N3)*TEMP 


D 330 


0296 




DO 119 K=l,3 


0 331 


0297 




GO TO (U2fll3fll3)f K 


D 332 


0298 


112 


L=0 


D 333 


0299 




GO TO HA QQ 
L=NINT(!;-U-1 0C7 


0 334 


0300 


113 


D 335 


0301 


114 


H=NINT(K)-L 


0 336 


0302 




R=RR(KI 


0 337 


0303 




DO 115 N=lfH 


D 338 


0304 




ITEMPflNiL 


0 339 


0305 




GX«N)*5UmiTEMP)*R 


0 340 


0306 


115 


R=R+H|M 


D 341 


0307 




FP=FP+S|HP(H,GX,MiK) 


0 342 


0308 




IF (K.eQ.3) FP=FPfSUM(N3)*RR(4»«H(4)/3. 


D 343 


0309 




GO TO ( lUfU7,ll8)t K 


0 344 


0310 


116 


FPR(K)'sPP 


D 345 


0311 




CO TO U9 


D 346 


0312 


117 


fPRIKJ*|FP-FPR(K-l)) 


.. . J., 0 347 
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0 348 
0 349 
0 350 
0 351 
0 352 
0 353 
0 354 
0 355 
D 356 
0 357 
0 358 
0 359 
0 360 
D 361 
0 362 
0 363 
0 364 
0 365 
n 366 
0 367 
D 368 
0 369 
0 37C 
D 371 
D 372 
D 373 
D 374 
0 375 
D 376 
D 377 
D 378 
D 379 
0 380 
0 381 
D 382 
n 383 
0 384 
0 385 

C • ■ 0 386 

C CALCULATION OF POWER NORMALIZED KUXES J ^JJ 

D 38B 

DO 131 1=1.3 I III 

00 131 N=1,N4 I 
FLUX(I,N) = FLOX !! «i5AM ^ JJJ 

AlDD{IrN)=FLLXlM' ■ " f:: 



0313 




GO TIj 119 






1 1 o 


FPR{K)=(FP-FPRIK-l)-FPR<K-2)) 








CQf^TINUE 








FP^OoO 








DO 120 K=li3 




1 ft 




FPR(K)=2#*P1*FPI<(K) 






120 


fpzFP-fFPRIK) 








GAM=POW/FP 








T£KP=PI/SQRT(BH2) 








GAK=6AM/TEMP 








WRITE i6il62) 








WRITE (6vl63) (RR( L) iL=lt3) 








WRITE (6fl64) EIGEN 








WRITE (6)165) 








GO TO (121f 121f 121f 122il23,124)i 


ICP 


\JJC0 


121 


WRITE (61I66) ICP 








GO TO 125 






122 

4 fa€* 


WRITE (6«167) 




C33i 




GC TO 125 






1 23 

X b 


WRITE I6il68) 








ruj TO 125 




U J J" 


124 


WRITE (6fl69) 






125 


DO 126 K=lf3 




0336 


126 


SCA(3fK>=SCA(3iK)-DI3tK)*BH2 




V J J f 




GO TO (127,128il29il27f I27il28) f 


ICP 




127 


K = l 








GO TC 130 




P^4'' 


128 


K=2 








GO TC 130 




r.A4? 


129 


K=3 




n^4A 


1 30 


TEMP'SCAI 3,K)-SCA0( 3»K) 








WRITE (6tl70) TEHP 




U3h!> 




TFMPsTFMP /SCP*1 .£+24 








toRfTf (6il71) TEMP 




0347 




WRITE (6,172) 




0348 




WRITE (6,173) {FLUX(l,N),N=l.N4) 




0349 




WRITE (6,174) (FLUX(2.N),N=1,N4) 




C350 




WRITE (6,175) (FLUX(3,N),N=1,N4) 





C 



0351 
C352 
0353 

0354 131 



0355 WRITE (6,176) 



0356 
0357 
0358 

0359 DO 132 Ll=l,3 

0363 LP=LI 



0361 

0362 132 



D 393 



WRITE (6,177) (AI0D(1«N),N=1,N4) D 
WRITE (6,178) (AI0D(2,N),N=1,N4) q\\ ^ 

WRITE (6,179) (AI00(3,N),N=1,N4) " ^ 



ERIC 



CALL PLOT (LI,N4,LP) ^ 
CONTINUE X 

0363 WRITE (6,130) 

0364 DO 133 1=1,3 

0365 TEMP=FLUX(1,1) 

0366 00 133 N=1,N4 

0367 133 AI00(I,N»=FLUX1I,N)/TEMP 
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WRITE (6,181) (AI0D(1,N),N=1,N4) J 406 

WRITE (6,182) (AIDD(2,N),N=1,N4I 0 ^07 

WRITE (6,183) (AI0D(3,N),N=1,N4) D 

DO 134 Ll=l,3 I 

LP=LIf3 ° 

CALL PLOT (LI,N4,LP) ^ 

134 CONTINUE J 
00 135 1=1,2 0 l\l 
DO 135 N=1,N4 J )JJ 

135 AIDD(I,N)=FLIJX(I,N)/FLUX(3,N) J 415 
WRITE (6,184) ° Jj; 
WRITE (6,185) (AIDD(1,N),N=1,N4) D ^IJ 
WRITE (6,186) (AI0D(2,N),N=1,N4) 0 2,o 

DC 136 Li=i,2 ° 7i; 

LP^^HI 2 4.0 

CALL PLOT (LI,N4,LP) " 2,, 

136 CONTINUE [ ill 
00 137 N=l,N4 0 423 

137 AIDD{l,N)=FLLX(l,N)/FLUX(2,N) D 424 
WRITE (6,187) (AIDD(1,N),N=1,N4) 0 425 
CALL PLOT (1,N4,9) g J26 

C CALCULATION OF POwER GENERATED IN EACH REGION 0 428 

C ^ 

DO 138 K=l,3 I J^O 

TEMP=PI/S0RT(eH2) ^ J^J 

FPR(K)=FPR(K)*TEMP*GAH ^ 2 ^ 

CONTINUE ^ y.l 

WRITE (6,188) (FPMK),K=1,3) ^ 
00 139 K=l,3 ill 

FPR(K)=FPR(K)/TEMP J" 

TEMP=PI»(RR(K+l)*PR(K + l)-RR(K)*RR(K)) " '»3/ 

FPfi(Kl=FPR(K)/TEHP ^ J38 

0398 139 CONTINUE ^ "ri 

WRITE (6,139) (FPR(K),K=1,3) " J^? 

WRITE (6,190) CVREL " 

IF (ICC.GT.O) GC TO 1 ^ 



0368 
0369 
0370 
0371 
0372 
0373 
0374 
0375 
0376 
0377 
0378 
0379 
0380 
0381 
0382 
0383 
0384 
0385 
0386 
0387 
0383 



0389 
0390 
1..91 

039'' 138 
0j93 
0394 
0395 
C396 
0397 



0399 
0400 
04J1 

C402 GO TO 141 



D 443 

0403 140 WRITE (6;i91) ^ 

0404 141 STOP nil 



0406 143 

0407 144 

0408 145 

0409 146 

0410 147 



ERIC 



0405 142 FORMAT (6EI0.0,6l 1X,12) ) Al S 111 

FORMAT (5E10.0) J III 

FORMAT (3E10.0) ^ 

FORMAT (4E10.0) ^ 

FORMAT (F5.3) " 

uti^ i-r. FORMAT (IHl) , , n fit 

0411 148 FORMAT (16H PROGRAM WAN DIC,///,26H REGIONS IN OP"-R 1,2,3,4,//) 0 453 

0412 149 FORMAT (17H RADII OF REGIONS, //,2X,8HREGI0N 1,4X 8HREG10N 2,4X,8HR D 454 

lEGION 3,4X,8HREGI0N 4,4X,8HREGI0N 5,//,2X,5(F8.2,4X)) 0 455 

0413 150 FORMAT (7h NINT =,4( 12, IX) , /,7H 0VREL=,F6.2,6X,6H PREC=,E9.2,6X,6ri 0 456 

0414 151 format' (//,6H BH2 =,E12.5,6X,28H POISON MICRO CR0SS-SECT10N=,IX,El 0 458 

12 • 5 / /) ^ ^ 

0415 152 FORMAT (26H ABSORPTION CROSS-SECTIONS, //,5X,7HGROUP l,7X,7HGRQUP 2 0 460 

l,7X,7HGR0Uf 3,/,3(2X,E12.5),/) ^ 

0416 153 FORMAT (//,23H DIFFUSION COEFFICIENTS, //,3(2X,E12. 5) ,/) 0 462 

0417 154 FORMAT (//,26H TRANSFER CROSS-SECTIONS ,//,3(2X,El2.5) ,/ 1 D 463 
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MAIN 



DATE = 74327 



12/56/19 



0413 


155 


FORMAT 


(//,23H NU-FISSION ,//i3(2X,E12.5I ,/) 






0419 


156 


FORMAT 


(//il4H KAPPA-FISSION>//f 342X,E12.5),/) 






C4?0 


157 


FORMAT 


(13H POWER (W) = ,E12.5) 








158 


FORMAT 


(IHL) 






0422 


159 


FORMAT 


(lXt4H K=,F10.6) 








160 


FORMAT 


(//,5H XP= ,E12.5,5X,15H EIGENVALUE K =,Ei2.5) 






04P4 


161 


FORMAT. (7H ITERs ,12) 






0425 


162 


FORMAT 


(Ihl.BH RESULTS,///) 






0426 


163 




I17H RADII OF REGIONS, //,2X,8HRtGI0N 1,4X, SHREGIbN 


2 


i4X,eHR 






,4X,8HREGI0N 4,4X,8HREGI0N 5.//,2X,5(F8.2,4Xl ) 






0427 


164 


FORMAT 


(//,4lh EIGENVALUE OF THE UNPOISONED REACTOR K =,E12.5) 




165 


FQRHAT 


(//,49H CONTROL POISON HAS USED IN THE FOLLOWING REGIONS,/) 


0429 


166 


FdRMAT 


(9H REGinN tI2) 






0430 


167 


FCRHAT 


i?\H All FISSION REGIONS ) 








168 




i 17H REGIONS 1 AND 2 ) 






0432 


169 


FORMAT 


( 17H RFGIQNS 2 AND 3 ) 






0433 


170 


format 


(//•36H MACRO CROSS^SECT ION OF THE POISON= tE12.5l 






0434 


171 


FORMAT 


(//*74H POISON CONCENTRATION = •E12.5) 






0435 


172 


PriDMAT 
r uwnft 1 


(////, 28H FLUXES NORMALIZED TO SOURCE) 






0436 


17S 


FORMAT 


(//,19H FAST FLUX-GROUP 1 ,/,5(2X,E12.5l ) 






0437 


174 


FORMAT 


|//,24H RESONANCE FLUX-GROUP 2 ,/,5(2X,E12.5) > 






043B 


175 


FORMAT 


(//,19H SLOW FLUX-GROUP 3 » / , 5( 2X,E12.5) ) 






C439 


176 


pnRMA T 

ru»^nfl 1 


(////, 28H FLUXES NORMALIZED TO PCHER ) 






0440 


177 


cnjM AT 


(//,19H FAST FLUX-GROUP 1 ,/,5{2X,El2.5) ) 






0441 


178 


? OR MAT 


(//,24H RESONANCE FLUX-GROUP 2 , /, 5(2X, £12.5) ) 






044P 


179 


FORHAT 


(//,19H SLOW FLUX-GROUP 3 ,/,5(2X,E12.5) ) 






0443 


180 


. FCRMAT 


t///,l6H RELATIVE FLUXES) 






0444 


181 


FORMAT 


(//,?H GROUP 1 ,/.5(2X,tl2.5)) 








182 


FORMAT 


(//,9H GROUP 2 ,/,5(2X,E12.5)) 








183 


FORMAT 


(//,9H GROUP 3 ,/,5(2X,E12.5)) 








184 


FORHAT 


(///,35H RATIOS OF FLUXES TO THE SLCW FLUX ) 








185 


FORMAT 


(//,19h GROUP 1 / GROUP 3 ,/,5(2X,£l2.5) ) 








I QO 


FORMAT 


(//»19H GROUP 2 / GROUP 3 ,/, 5(2X, E12. 5) ) 








187 


FORMAT 


(//,19H GROUP 1 / GPOUP 2 ,//,5(2X,E12.5)) 






045L 


188 


FORMAT 


(//i3lH POWER OUTPUT IN REGIONS I,2i3,//,2X,3(E12.5 


,5X)) 


0452 


189 


FORMAT 


(//,42H POKER DENSITY IN REGIONS 1,2,3 (K^/LITER), 


// 


,2X,3(E 






112. 5, 5X)) 






0453 


190 


FORMAT 


l//,33H COMPUTED OVEPRELAXATION FACTOR^ ,F6.3) 






0454 


191 


FORMAT 


(22H HI DF ITERATIONS > 400) 






0455 




END 









D 464 
D 465 
0 466 
0 467 
468 
469 
470 
471 
472 
D 473 
D 474 
475 
476 
477 
478 
479 
480 
D 481 
0 482 
0 483 
484 
485 
4B6 
487 

489 
0 490 
r 491 
[! 492 
493 
494 
D 495 

0 496 
D 497 
I) 498 

1 499 
I) 500 
') 501 
D 502 
D 503 



D 
D 
D 
D 
0 



D 
D 
D 
0 
0 
D 



D 
D 
D 
0 
•) 
P 



D 



42 



ERIC 



41 



PftOGRAH WAN QIC 

REGIONS IN ORDER Iv2f3f4f 

RAOI t OF REGIONS 

REGION 1 REGION 2 REGION 3 

0^0 87.70 124.00 

NINT =30 20 20 4 

0VREL= 1.90 PREC= O.lOE-03 



REGION 4 REGION 5 
152.00 159.00 
CONV== O.lOE-01 



6H2 = 0.73800E-04 POISON MICRO CROSS-SECTION^ 0.21090E 04 

ABSORPTION CROSS-SECTICNS 



GROUP 1 
0.35410E-02 

0.35410E-02 

0.35410E'02 

0*36220e-03 



GROUP 2 
0.25680E-01 

0.25680E-01 

0.25680E-01 

0.57220E-03 



GROUP 3 
0.13870E 00 

0.13870E 00 

6.13870E 00 

0.84730E-02 



DIFFUSION COEFFICIENTS 

0.17070E 01 0.85360E 00 0.39510E 00 

0.17070E 01 0.85360E 00 0.39510E 00 

0.1707DE 01 0.85360E 00 0.39510E 00 

0.19300E 01 0.91730E 00 C.31470E 00 

TRANSFER CROSS-SECTIONS 



0.61130E 


-01 


0.49440E- 


01 


0.0 


0.61130E 


-01 


0*49440E- 


01 


0.0 


0.61130E 


-01 


0.49440E- 


01 


C«0 


0.26280E 


-01 


0.28110E- 


01 


0.0 



NU-FISSIUN 
0.42920E-02 
0.42920E-02 
0.42920E-02 



0.15060E-01 
0.15060E-CI 
0.15060E-01 



0*16400E 00 
0.16400E 00 
0.1640GE 00 



EKLC 



KAPPA-FISSION 
0.93280E-13 
0.93280E-13 
0.93280E-13 

POMER (Wl 



0.19010E-12 
0.19010E-12 
0.19010E-12 
0.24410E 10 



0.21290E-11 
0.21290E-11 
0.21290E-I1 
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RESULTS 



RAOII OF REGIONS 



REGION 1 



REGION 2 



REGION 3 



REGION 4 



REGION 5 



0.0 



87.70 



152.00 



159.00 



EIGENVALUE OF THE UNPOISONED REACTOR K = 0.97971E 00 

CONTROL POISON MAS USED IN THE FOLLOWING REGIONS 
ALL FISSION REGIONS 

MACRO CROSS-SECTION OF THE POlSCN= -0.37987E-02 
POISON CONCENTRATION = -0.16012E 19 



FLUXES NORMALIZED TO SOURCE 



FAST FLUX-GROUP 1 



0.^2922E 


-03 


0.4290<»E- 


-03 


0.42848E- 


-03 


0.42755E 


-03 


0.42625E 


-03 


0.^2^57E 


-03 


0.^2253E 


-03 


0.42012E- 


-03 


0.41734E 


-03 


0.41420E 


-03 


0.^1070E 


-03 


0.40683E- 


-03 


C.40261E 


-03 


0.39804E- 


-03 


0.39312E 


-0 3 


0.38785E 


-03 


0.3822^E 


-03 


C.37629E 


-03 


0.37002E 


-03 


0.36 342E' 


-03 


0.35650E 


-03 


0.3^927E 


-03 


0.34173E 


-03 


0.33390E 


-03 


0.32578E 


-03 


0.31737E 


-03 


0.30870E- 


-03 


0.29977E 


-03 


0.29059E 


-03 


0.28117E 


-03 


0.271535 


-03 


0.2654 3E 


-03 


0.25925E 


-03 


0.25300E 


-03 


0.24667E 


-03 


0.2^327E 


-03 


0.23381E-03 


0.;i2728E- 


-03 


0.22069E 


-03 


0.21404E 


-03 


0.2073^E 


-03 


0.20059E 


-03 


0.19379E 


-03 


0.18696E 


-03 


0.18008E 


-0 3 


0.17317E 


-03 


0.L6623E 


-03 


C.15926E 


-03 


0.15228E 


-03 


0.14527E 


-03 


0.13825E 


-03 


0.132d2E 


-03 


C.12740E 


-03 


0.12197E 


-03 


0.11654E 


-03 


O.lllUE 


-03 


0.10569E 


-03 


C.10027E 


-03 


0.94856E 


-04 


0.d9454E 


-04 


0.8^063£ 


-0^ 


0* 78686E 


-04 


0.73325E 


-04 


0.67980E 


-04 


0.62653E 


-04 


0.573^1E 


-0^ 


0.52038E 


-04 


C.46722E 


-04 


0.41343E 


-04 


0.35784E 


-04 


0.2977^E 


-0^ 


0.23213E 


-04 


0.17707E 


-04 


0.13012E 


-04 


0.89179E 


-05 



RESONANCE FLUX-GROUP 2 



0.i4819E 


-03 


0.34804E- 


-03 


C.34759E 


-03 


0.34683E 


-03 


0.34577E 


-C3 


0.34441E 


-03 


0.34275E 


-03 


C.34080E 


-03 


0.33854E 


-03 


0.33 599E 


-03 


0.33314E 


-03 


0.33001E 


-03 


0.32658E 


-03 


0.32287E 


-03 


0.31887r 


-0 3 


0.31459E 


-03 


0.31004E 


-03 


G. 3052 IE 


-03 


0.30012E' 


-03 


0.29476E 


-03 


0.28914E 


-03 


0.2832 7E 


-03 


0.27716E 


-03 


0.27080E 


-03 


0.26421E 


-03 


0.25739E 


-03 


U.25036E 


-03 


0.24311E 


-03 


0.23566E 


-03 


0.22802E 


-03 


0.22019E 


-03 


0.21525E 


-03 


0.2 102 3E 


-03 


0.20516E 


-03 


0.23003E 


-0 3 


0.19484E 


-03 


0.18959E 


-03 


0.18429E 


-03 


0.17895E 


-03 


C.17356E 


-03 


0.16812E 


-03 


0.16265E 


-03 


0.15714E 


-03 


0.15159E 


-03 


0.14602E 


-03 


0.14041E 


-03 


0.13478E 


-03 


0.12913E 


-03 


0.12347E 


-03 


0.11778E 


-03 


0. 11209E 


-03 


0.10769E 


-03 


0.10329E 


-03 


0.98890E 


-04 


0.94487E 


-04 


0.90086E 


-04 


0.85687E 


-04 


0.81293E 


-04 


0.76904E 


-04 


0.72522E 


-04 


0.6ei49E 


-04 


0.63786E 


-04 


0.5943iE 


-04 


0.55091E 


-04 


0.50761E 


-04 


0.46441E 


-04 


0.42132E 


-04 


C.37833E 


-04 


0.33548E 


-04 


0.29291E 


-04 


0.25108E 


-04 


0.20280E 


-04 


0.15418E 


-04 


0. 105i7E 


-04 


0.55S11E 


-05 



SLOW FLUX- 


GROUP 


3 
















0.1274BE 


-03 


0.12742E 


-03 


C.12726E 


-03 


0.12698E 


-03 


0.12659E 


-03 


0.12609E 


-03 


0.12549E 


-03 


C.12477E 


-03 


0.12394E 


-03 


0. 12 3C1E 


-C3 


0.12197E 


-03 


0.12082E 


-03 


0.11956E 


-03 


0.1182UE 


-03 


0.11674E 


-03 


0.11517E 


-03 


0.11351E 


-03 


C.11174E 


-03 


0.1G987E 


-03 


0.10791E 


-03 


0.10586E 


-03 


0. 10371E 


-C3 


0.10147E 


-03 


0.99139E 


-04 


0.96725E 


-04 


0.94229E 


-04 


0.91653E 


-04 


0.89000E 


-04 


0.d6272E 


-04 


C.83474E 


-04 


0.806C9E 


-04 


0.78797E 


-04 


0.76962E 


-04 


C.75105E 


-04 


0.73225E 


-04 


0.71325E 


-04 


0.69405E 


-04 


C.67466E 


-04 


0.65509E 


-04 


0.63535E 


-04 


0.61545c 


-04 


0.59541E 


-04 


C.57524E 


-04 • 


0.55494E 


-04 


0.53452E 


-C4 


0.51401E 


-04 


0.49340E 


-04 


0.47272E 


-04 


0.45197E 


-04 


0.43117E 


-04 


0.41033E 


-04 


0.39423E 


-04 


0.37812E 


-04 


0.36200E 


-04 


0.34588E 


-04 


0.32977E 


-04 


0.31367E-0^ 


0.29758E 


-04 


0.28152E 


-04 


0.26548H 


-04 


0.24948E 


-04 


0.23353E 


-04 


0.21764E 


-04 


0«20186E 


-04 


0.18625E 


-04 


0. 17100E 


-04 


0.15647E-04 


0.14350E 


-04 


0.13391E 


-04 


0.13176E 


-04 


0.14593E 


-04 


0.16797E 


-04 


0. 14841E-04 


0.99331E 


-05 


0.30294E 


-05 



ERIC 
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FLUXES NORMALIZED TO POWER 



FAST FLUX-GROUP 
0.21237E 15 
0.21007E 15 
0.20320E 15 
0.19190E 15 
0.17639E 15 
0.15703E 15 
0.13435E 15 
0.11888E 15 
0.10259E 15 
0.85681E 14 
0.68402E 14 
0.54975E 14 
0.41593E 14 
0.28371E 14 
0.14731E 14 



1 

0.21228E 15 
0.20906E 15 
0.20129E 15 
0.18912E 15 
0.17281E 15 
0.15274E 15 
0.13133E 15 
0.11568E 15 
0.99247E 14 
0.82248E 14 
0.65719E 14 
0.52291E 14 
0.38932E 14 
0.25747E 14 
0.11485E 14 



0.21200E 15 
0.20787E 15 
0.19920E 15 
0.18618E 15 
0.16908E 15 
0.14832E 15 
0.12827E 15 
0.11245E 15 
0.95885E 14 
0.78801E 14 
0.63034E 14 
0.49610E 14 
0.36279E 14 
0.23117E 14 
0.87612E 13 



0.21154E 15 
0.20649E 15 
0.19694E 15 
0.18308E 15 
0.16520E 15 
0.14378E 15 
0.12518E 15 
0.10919E 15 
0.92502E 14 
0.75343E 14 
-0.60347E 14 
0.46933E 14 
0.33635E 14 
0.20456E 14 
0.64379E 13 



0.21090E 15 
0.20494E 15 
0.19451E 15 
0.17981E 15 
0.16119E 15 
0.13912E 15 
0.12205E 15 
0.10590E 15 
0.89100E 14 
0.71876E 14 
0.57661E 14 
0.44260E 14 
0.30999E 14 
0.17705E 14 
0.44124E 13 



RESONANCE FLUX-GROUP 2 



0. 


17228E 


15 


0.17220E 


15 


0.17198E 


15 


0.17160E 


15 


0.17108E 


15 


0. 


17041E 


15 


0.16959E 


15 


0.16862E 


15 


0.16750E 


15 


0.16624E 


15 


0. 


16483E 


15 


0.16328E 


15 


0.16159E 


15 


0.15975E 


15 


0.15777E 


15 


0. 


15565E 


15 


0.15340E 


15 


0.I5101E 


15 


0.14849E 


15 


0.14584E 


15 


0. 


14306E 


15 


0,14016E 


15 


0.13713E 


15 


0.13399E 


15 


0.13073E 


15 


0. 


12735E 


15 


0.12387E 


15 


0.12029E 


15 


0.11660E 


15 


0.11282E 


15 


0. 


10895E 


15 


0.10650E 


15 


0.10402E 


15 


0.10151E 


15 


0.98969E 


14 


0. 


96400E 


14 


0.93805E 


14 


0.91185E 


14 


0.88540E 


14 


0.85872E 


14 


0. 


83Id4E 


14 


0.80475E 


14 


0.77748E 


14 


0.75005E 


14 


0.72246E 


14 


0. 


69473E 


14 


0.66688E 


14 


0.63893E 


14 


0.61089E 


14 


0.58277E 


14 


0. 


SS-^bOE 


14 


0.53284E 


14 


0.51107E 


14 


0.48929E 


14 


0.46750E 


14 


0. 


445 72 E 


14 


0.42396E 


14 


0.40222E 


14 


0.38050E 


14 


0.35882E 


14 


0. 


33719E 


14 


0.31560E 


14 


0.29406E 


14 


0.27258E 


14 


0.25115E 


14 


0. 


22978E 


14 


0.20846E 


14 


0.18719E 


14 


0.16599E 


14 


0.144'>3E 


14 


0.I2423E 


14 


0.10034E 


14 


0.76287E 


13 


0.52137E 


13 


0«276. E 


13 



SLOW FLUX-GROUP 3 



0.63073E 


14 


0.63046E 


14 


0.62963E 


14 


0.62826E 


14 


0.62635E 


14 


0.62388E 


14 


0.62088E 


14 


0.61733E 


14 


0.61324E 


14 


0.60862E 


14 


0.60347E 


14 


0.59778E 


14 


0.59158E 


14 


0.58485E 


14 


0.57761E 


14 


0.56986E 


14 


0.56161E 


14 


0.55286E 


14 


0.54363E 


14 


0.53392E 


14 


0.52375E 


14 


0.51311E 


14 


0.50203E 


14 


0.49052E 


14 


0.47858E 


14 


0.46623E 


14 


0.45348E 


14 


0.44035E 


14 


0.42686E 


14 


0.41301E 


14 


0.39883E 


14 


0.38987E 


14 


0.38079E 


14 


0.37160E 


14 


0.36230E 


14 


0.35290E 


14 


0.34340E 


14 


0.33380E 


14 


0.32412E 


14 


0.31436E 


14 


0.30451E 


14 


0.29460E 


14 


0.28461E 


14 


0.27457E 


14 


0.26447E 


14 


0.25432E 


14 


0.24412E 


14 


0.23389E 


14 


0.22 363E 


14 


0.21333E 


14 


0.20302E 


14 


0.19506E 


14 


0.18708E 


14 


0.17911E 


14 


0.17114E 


14 


0. 16316E 


14 


0.15520E 


14 


0.14724E 


14 


0.13929E 


14 


0.13135E 


14 


0.12344E 


14 


0.11554E 


14 


0.10768E 


14 


0.99874E 


13 


0.92153E 


13 


0.84606E 


13 


0.77418E 


13 


0.70999E 


13 


0.66258E 


13 


0.65193E 


13 


0.72204E 


13 


0.83110E 


13 


0.73431E 


13 


0.49147E 


13 


0. 14989 E 


13 



ERIC 
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PLOT OF FAST FLUX - GHOUP I 



0.^60E 1<> 



Z-AXIS 
0.676b 1^ 



0.1246 IS 



0.171F l& 



C.?i2f 15 



<-AXIS 
0.0 

0.29^E 01 

0.586E 01 

0.e77£ 01 

0.U7E 0^ 

0.U6E 02 

0.l7bE 0^ 

0.23^E 0^ 

0.2oi£ 02 

0.29^E 02 

0«322(: 02 

0.3:>1E 02 

0.380E 02 

U.409C 02 

0.4ia£ 02 

0.4<»7E 02 

0.t>2(>E Ji! 

0.555= 02 

0.5d5E 02 

O.bl^F 02 

u.6^jt 02 

0.6726 02 

0.702E 02 

0.73iE 02 

0.7O0E 02 

0,7696 02 

0.dl9E 02 

0.a^3E 02 

0.87/'^ 02 

U.a95E 02 

0.9liE 02 

0.931E 02 

0.95ut: 02 

0.96ar 02 

0.9JO? 02 

U.LOjF 03 

0.1j2c J3 

O.KHE 03 

0.106E 03 

O.lOtiF 03 

0.1U9E 05 

O.lllE 03 

0.113E 03 

0.115E 03 

0.117E 03 

0.119E u3 

0.120E 05 

o•l^^f Oi 

0.12^(- 03 

0.125E 33 

U.12/E 03 

0.1283 03 

0.130E 0 3 

U.131E 03 

0.132E 03 

0.l3<>c 03 

0.135E 03 

0.137H 03 

0.138E C3 

0,139E 03 

0.141H 03 

0.U2C 03 

O.U<>E 03 

u.Ut>E 03 

O.K6t 03 

0.U8E 03 

0.U9F 03 

0.151E 03 

0.152E 03 

O.li^E 33 

0.155E C3 

0.157E 03 

0.159E C3 



.......... 



ERIC 
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PIOT OF RtSONANCE FLUX -GROUP 2 



2-AXlS 

^•^^*>^ 13 0.367E 0.706t O.iO^e 15 0.138C lb L.l72c 15 

K^AXIS 

o»a » 

0.t»ti3E 01 » 
0.677E 01 « 
0.117E 02 « 
0.i^6E 02 * 
0.175E 02 » 
0.20*>E 02 <y 
0.23Ae 02 IT 
0.203E 02 « 
0*292? 02 * 
0.3226 02 * 
0.3516 02 , 
0.3dCE 02 « 
0.^09€ 02 « 
U.43ti£ 02 , 
O.^bdc 02 * 
0.<i<^rE 02 « 
0.526e 02 
0.555E 02 

0.5d5E 02 * 
O.ol^E 02 « 
0.6^iE 02 * 
J.6I2E 02 . 4 

0.702C 02 « 
0.731E 02 « 
0.760E 02 « 
0.789E 02 ♦ 
U.eiSE 02 p 
O.b^HE 02 * 
U.b77E u2 ^ 
0.695= 02 p 
0.913^ 02 * 
0.931J: 02 ♦ 
0.950'? 0^ « 
0.96eE 02 * 
O.Stl6E 02 * 
O.lJOE 0) « 
0.I02E 03 4. 
U.lO^E C3 « 
O.lOfiE 03 « 
O.lOaE 03 « 
0.109E 0 3 41 
O.lUe 03 • 
0.1l3e 03 m 
0.115E 03 « 
0.li7E 0:t * 
0.119E 03 » 
0.120E 03 « 
0.122E 03 * 
O.l^^E 03 i» 
J.125E 03 9 
U.127E 03 « 
0.128^ C3 * 
0.130E 03 * 
0.131F 03 * 
.0.132E 03 ♦ 
0.134E 03 

0.1J5E 03 * 
0.137E 03 ♦ 
0.138E 03 m 
-©•139E 03 * 
0.I41E 03 * 
0.142E 03 m 
0.144E 0 3 * 
0.I45E 03 ♦ 
0.146^ 03 • 
0.l4aE 03 • 
0.i49E 03 * 
0.151F 03 * 
0.152E 03 * 
0.154E 03 • * 

0.15*>E 03 * 
0.157E J3 » 
0.159E 03 ♦ 



EKLC 
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PLOT OF SLOW FLUX - GKUUP 3 



0.15UE 13 



0.138E 14 



Z-AXIS 
0.26U 14 



0.384E 14 



0.50HE 14 



0.6311 14 



01 



02 
02 



02 
02 



0.2<^2E 

U.5d3R 01 

0.677E 01 

0.117F 02 

0.146E 02 

0.175E 02 

0.205E 02 

0.234E 02 

0.2&3E 02 

0.292E 02 

0.322E 02 
0.351E 
0.3805 

0.40Vt: u2 

0.43Be U2 
0.46d'' 
0.497? 
U.526E 

0.5^)56 j2 

0.385E 02 

0.614E 02 

0.643E C2 

U.672E 02 

O.702E 02 

0.731E 0^ 

J.766E 02 

0.7d9E 02 

0.019^ 0^ 

O.840E 02 

0.b77E 02 

0.89t>E 02 

0.91 ^E 02 

0.93lt 02 

0.9S0e 02 

0.96tiF 02 

O.S06E 02 

O.lOOt 03 

0.102E 03 

0.104c 03 

0.106E -13 

O.lJaE 03 

0.1 J9E Oi 

O.lllE 03 

0.113E 03 

0.115F 03 

0.117<^ 03 

j.ll9E 03 

0.120E 03 

U.122E 03 

0.1 24t 0 3 

0.125E 03 

0.127E 03 

0.126E 03 

0.130E 03 

0.131E 03 

0.1326 03 

0,134E 03 

0.135E 03 

0.137E C.V 

a.l3dE: 03 

0.U9E 03 

0.141E 03 

0.142E 03 

0.144E 03 

0.145E 03 

0.l4bE 03 

0.146E 03 

0.U9E 03 

0.151E o:) 

0.152c 03 

0.154E 03 

3.155E 0^ 

0«137E 03 

0«159E 03 



EKLC 



8 



RELATIVE FLUXES 



- 47 - 



GROUP 1 



O.IOOOOE 


01 


0.99957E 


00 


0.99826E 


00 


0.99609E 


00 


0.99306E 


00 


0.98916E 


00 


0.9B440E 


00 


0.97879E 00^ 


fOi9725?E/00 
0.9?735E^00 


0.96500E 


00 


0.95683E 


00 


0.94783E 


00 


0.93800E 


00 


0.91588E 


00 


0.90360E 


00 


0.89053E 


00 


0.87668E 


00 


0.86206E 


00 ' 


0. 84668 E 


00 


0.83056E 


00 


0.81371E 


00 


0.79616E 


00 


0.77791E 


00 


0.75899E 


00 


0.73942E 


00 


0.71921E 


00 


0.69841E 


00 


0.67702E 


00 


0«65507E 


00 


0.63260E 


00 


0.61839E 


00 


0.60400E 


00 


0.58943E 


00 


0.57469E 


00 


0.55978E 


00 


0.54472E 


00 


0.52950E 


00 


0.51415E 


00 


0.49 866E 


00 


0.48305E 


GO 


0.46733E 


00 


0.45 150E 


00 


0.43557E 


00 


0.41955E 


00 


0.40345E 


00 


0.38728E 


00 


0.37105E 


00 


0.3M77E 


00 


0.33 845E 


00 


0.32209E 


00 


0.30945E 


00 


0.2968 IE 


00 


0.2a416E 


00 


0.27151E 


00 


0.25886E 


00 


0.24623E 


00 


0.23360E. 00 


0.22099E 


00 


0.20841E 


00 


0.19585E 


00 


0.18332E 


00 


C.17083E 


00 > 


0.15838E 


00 


0.14597E 


00 


0.13359E 


00 


0.12124E 


00 


0.10885E 


00 <; 


P.96321E-01 


0.83369E- 


-01 


0.69366E' 


-01 


0.54081E- 


-01 


0.41254E-01 


' 0.30314E- 


-01 


0.20777E- 


-01 



GROUP 2 

O.IOOOOE 01 
0.98915E 00 
0.95679E 00 
0.90351E 00 
0.83042E 00 
0.73923E 00 
0.63239E 00 
0.55957E 00 
0.4a285E 00 
0.40326£ 00 
0.32192E 00 
0. 258736 00 
0.19572E 00 
0.13338E 00 
0.72109E-01 



0.99957E 00 
0.98439E 00 
0.94778E 00 
0.89043E 00 
0.81356E 00 
0.71902E 00 
0.61818E 00 
0.54450E 00 
0.46713E 00 
0.38710E 00 
0.30930E 00 
0.24609E 00 
0.18319E 00 
0.12100E 00 
0.58245E-01 



0.99826E 
0.97876E 
0.93794E 
0.87657E 
0.79600E 
0.69821E 
0.60379E 
0.52929E 
C.45130E 
0.3708 7E 
0.29666E 
0.2334 IE 
0.17069E 
0.10866E 
0.44281E- 



00 0.99609E 00 
00 0.97229E 00 
00 0.92727E 00 
00 0.86194E :&0 
00 0^7Z7.4£:^ 
00 0.6768^E;0b 
00 ,.0-58,9>22E 00 
00 .;;?lD.5i|94E 00 

00 4^0. 43^3 7E 00 

00 1 ' 0.35460E 00 
00 •0*28401 E 00 
00 0.22087E 00 
00 0.15822E 00 
00 0.96349E-01 
•01 0-30263E-01 



0.99305E 00 
0.96496E 00 
0.91580E 00 
0.84655E 00 
0.75881E 00 
0.65487E 00 
0.57447E 00 
0.49846E 00 
0.41936E 00 
0.33828E 00 
0.27137E 00 
0.20828E 00 
0.14578E 00 
0.84124E-01 
0.16029E-01 



GROUP 3 

O.IOOOOE 01 
0.98915E 00 
0.95678E 00 
0.90349E 00 
0.83038E 00 
0.73918E 00 
0.63234E 00 
0.55951E 00 
0.48279E 00 
0.40321E 00 
0.32188E 00 
0.25869E 00 
0. 19571E 00 
0.i3414E 00 
0.11448E 00 



0.99957E 00 
0.98438E 00 
0.94776E 00 
0.89041E 00 
0.81352E 00 
0.71897E 00 
0.61813E 00 
0.54445E 00 
0.46707E 00 
0.38705E 00 
0.30925E 00 
0.24606E 00 
0.18319E 00 
0.12274E 00 
0.13177E 00 



0.99826E 00 
0.97876E 00^ 
0.93792E 01,0 V 
0.87654E 00 
0.79596e 00 
0.69816E 00 
0.60373E 00 ' 
0.52924E 00 
0.45125E 00 
0.37083E 00 
0.29662E 00 
0.23344E 00 
0.17073E 00 
0.11257E 00 
0.11642E 00 



0.99'fe09E 00 
0.97228E 00 
0.92725E 00 
0.86191E 00 
^).77770E 00 
0.67677E 00 
0.58916E 00 
0.51388E 00 
0.43532E 00 
0.35455E 00 
0.28397E 00 
0.22084E 00 
0.15835E 00 
0.10505E 00 
0.77920E-01 



0.99305E 00 

0.96495E 00 

0.91577E 00 

0.84652E 00 

0.75877E 00 

0.65482E 00 

0.57442E 00 

0.49840E 00 

0.41931E 00 

0.33823E 00 

0.27133E 00 

0.20826E 00 

0.14611E 00 

0.10336E 00 
0.23764E-01 
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PLOT OF GROUP I 



3.20BC-01 



0.217E CO 



2-AXlS 
0.412E 00 



0.60dE 00 



0.804E 00 



O.irct CI 



0.0 

0.2*92c 
0««>65E 
0.877E 
0.117E 
0. 1460 
0. I75c 
0.2v}5e 

0.2631: 
0.2921: 
0.322P 
0.3J>1 = 

0.43tfF 

o*46be 

• W.4S/E 
0.52o6 
0.t)5*>b 

C.643E 
u.e72E 
:«.70,?z 
U.731: 

0. 700": 

0.7dSe 
0.H19F 

J.677- 

0.931r 

O.^doF 
u . I uCE 
0. iOZt 

a.ij4£ 

u. lOoE 
u. lOdE 
J.IOSE 
O.lllC 
J. U3E 
0. llbE 
0.117E 
0. llVr 
0. 120^ 
J.l22e 
0.124E 
J.125F 
0.i27r 
0.l2dr= 
0.130E 
0.131E 

. 0.132^ 
0.134E 
0. 1^5E 
0.i37b 
0.l3ttE 

' 0.13VE 
0.141E 
0.142E 
0. UnF 
0.U5C 
0.146L 

0. use 

0. U9E 
0.151E 
0.i52E 
3*1:>4E 
0.155E 
0.157s 
0.15SP 



01 
01 
01 
02 
02 
02 
02 
02 
02 
02 
02 
02 
02 
02 
02 
j2 
OJ 
1.2 
ji! 
02 
02 
02 
02 
02 
02 
j2 
02 
02 
OJ 

02 
u2 
02 
j2 
D2 
02 
03 
03 
03 
03 
03 
03 
03 
03 
03 
J3 
03 
33 
j3 
03 
03 
03 
■35 
03 
03 
03 
03 
03 
03 
03 
03 
03 
33 
03 
03 
03 
03 
C3 
03 
03 
?i 
03 
03 



o.i<tor 



J. Irt'J: 
v). 0.* 

.'.to-. »' 



J. ro^: 


0 : 


J. T iU 












J. Hl'*'^ 






OJ 




•).,' 




Jl* 


J'. 


J.' 


O.'lJil 










Jl* 








J 1 






j.u«.t 

u. UcL 


0 » 




J i 


. .lo>^^ 


^ 1 


J. tilt 


:i i 


j.i 1 j( 


OJ 


J.ll vt 


Ji 




ji 




Oi 


S..12JL 


u 1 




0 i 


J. w<>r 




J.1<!!>L 


Cll 


J.t.;Te 


0 i 


0.1 2a£ 


UJ 


O.ljOl 


3> 


O.lil' 


0 J 


0.1 i^l 


J} 


j.U*f 


u> 


u.l jVt 


0> 


j.i 


.1 i 


fc.l>HE 






J 1 



/-AXIS 
O.'tlQf 00 



3.l**s 

O.l^SE 

O.lh'il JJ 

J. 1^15 0» 

J.ltx.' 0\ 

j.lb1»£ 0) 

ii.lill JJ 

O.li^i: 3i 
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PLOT OF GROUP 2 

Z-AXIS 

0,l60b-0l 0.213E 00 O.^lOf 00 0.6066 00 0.803E 00 O.ICC^ CI 

• •*«•••••••••••••«••• 

R-AXIS 

0.0 * 

0.292C 01 * 

O.tibSE 01 * 

0.877E OL • 

0.1175 02 * 

0.U6E 02 * 

0.l75e 02 » 

0.20SE 02 *- 

0.234E 02 * 

0.263E 02 ♦ 

0.292i 02 * 

0.J2^5 02 • 

J.3bl5 C2 * 

0.380E 02 * 

O.^O^J? 02 * 

0.4336 02 • 

0.46dE 02 * 

0.497E J2 ♦ 

0.526t 02 * 

0.555.E 0^ * 

0.585c O'i * 

J.fcUE 02 * 

0.643^ 0^ * ' 

J.672E 0^ * 

0.702E 02 * 

J.73l£ 02 ♦ 

0.7oOE 02 * 

J.789E 02 • 

J.819E 02 * 

J. 8485 02 * 

U.d77P 02 » 

0.8955 0^ » 

0.9l^E 02 ^ 

0.931E D2 * 

0.9i0c 02 * 

0.«i68e 02 * 

0.9B6E 02 * 

J. 1005 OJ ♦ 

J.102E 0:i * 

J.1J4E 0 3 

U. ICcc 03 * 

u.lC8£ Oj ♦ 

C-.109E 03 * 

O.llLE 03 * 

0.U3P 03 * 

O.llbE 03 * 

0.U7E 03 • 

0.119E 03 • 

0.120E 03 * 

U.122E Oj • 

0.124E 03 * 

0.l25fc 03 ♦ 

0.127E 03 * 

0.128E OJ * 

0.130E 03 * 

O.ljlE 03 * 

0.132E 03 ♦ 
' 0.134E 03 * 

0.lJ5t Oj ♦ 

J. 137= 03 ♦ 

C.138E 03 • 
• 0.1J9E 0 3 * 

U.141E 03 * 

J.i42E 03 * 

0.144E Oj * 

0.14!>E 03 * 

0.146E 03 * 

0.148E (i3 ♦ 

0.149E 03 ♦ 

O.lblE 03 * 

0.152E 03 * 
03 * 

0.155E 03 * 

0.157E 03 * 

0.1!>9I: 03 • 
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PLOT OF GROUP 3 



t-AXIS 
U.O 

0.292E 01 
0.385E 01 
0.e77E 01 
O.IUP 02 
0.U6E 02 
U.175E 02 
0.205E 02 
0.234E 02 
0.263E 02 
0.292E 02 
0.322E 02 
0.351E C2 
0.380F 02 
0.<»OVE 02 
J.4JBE U2 
J.466E 02 
U. 4*^75 02 
0.526E 02 
O.SbSe 02 
0.58bE 02 
0.6UE 02 
0.643E 02 
0.672c 02 
0.7O2E 02 
0.731E 02 
U.760E v^2 
0.7U9E 02 
0.6L9E 02 
0.84bE 02 
0.677E 02 
0.d9bE 02 
0.913E j2 
0.93IG 02 
0.9505 02 
0.96BE 02 
0.*?S6E 02 
O.lJOE 03 
0.102F 03 
0.104E 03 
0.106^ 03 
O.lOtiE 03 
0.1096 Oj 
U.lllE 03 
O.I13!3 33 
U.115E Oi 
O.lUE 03 
O.llSE 03 
0.120E C3 
0.1225 Oj 
0.124E 03 
0.125E 03 
0.127": 03 
0.l2dE 03 
0.13GE 03 
0.131E 03 
0.132E 03 
0.134E 03 
W.135E 03 
0.137E 03 
0.13dE 03 
0.i39E 03 
O.UIE 03 
O.U2e 03 
0.U4E 03 
O.l^SE 03 
0.1^6E 03 
CUBE 03 
0.1^9E 03 
0.151E 03 
0.132E 03 
O.IS^E 03 
0.155E 03 
0.157F 03 
O.IS^E 03 



0.2idE-0l 0.219E 00 0.414E 00 0.610E 00 0.805E 00 O.IOCJ- 01 
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RATIOS OF FLUXES TO THE SLOW FLUX 



GROUP 1 / GROUP 3 



0* 33671E 


01 


0. 3367iE 


OX 


0.33671E 


01 


0.33671E 


01 


0.33671E 


01 


©•33671E 


01 


0.33671E 


01 


0.33672E 


01 


0.33672E 


01 


0.33672E 


01 


0.33673E 


01 


0.33673E 


OX 


0.33673E 


01 


0.33674E 


01 


0.33674E 


01 


0. 33675E 


01 


0.33675E 


OX 


0.33676E 


01 


0.33677E 


01 


0.33677E 


01 


0*336 /8E 


01 


0.33679E 


01 


0. 33679E 


01 


0.33680E 


01 


0.33680E 


01 


0.33&8IE 


01 


0.33682E 


01 


0.33683E 


01 


0.33683E 


01 


0.33684E 


01 


0.33685E 


01 


0.33685E 


01 


0.33686E 


01 


0.33686E 


01 


0.33686E 


01 


0.33687E 


01 


0.33687E 


01 


0.33688E 


01 


0.33688E 


01 


0.33688E 


01 


0.33689E 


01 


0.33689E 


01 


0.33689E 


01 


0.33690E 


01 


0.33690E 


01 


0.i3690E 


01 


0.33691E 


01 


0.3369 IE 


01 


0.. 33 69 IE 


01 


0.33692E 


01 


0.33692E 


01 


0.33692E 


01 


0.33693E 


01 


0.33693E 


01 


0.33693E 


01 


0.33693E 


01 


0.33694E 


01 


0.33694E 


01 


0.33694E 


01 


0.33695E 


01 


0.33695E 


01 


0.33695E 


01 


C. 3369 IE 


01 


0.33677E 


01 


0.33639E 


01 


0.33533E 


01 


0.33257E 


01 


0.32559E 


01 


0.30873E 


01 


0.27158E 


01 


0.20403E 


01 


0.13819E 


01 


0.11931E 


01 


0.13099E 


01 


0.29438E 


01 



GROUP 2 / GROUP 3 



0.27314E 


01 


0.27314E 


CI 


0.27314E 


01 


0.27314E 


01 


C.27314E 


01 


0.27314E 


01 


0.27314E 


01 


0.27314E 


01 


0.27314E 


01 


C.27314E 


01 


0.27314E 


01 


0.27314E 


01 


0.27314E 


01 


0.27314E 


01 


0.27315E 


01 


0.27315E 


01 


0.27315E 


01 


0.27315E 


01 


0.27315E 


01 


0.27315E 


01 


0.273156 


01 


0.27315E 


01 


0.27315E 


01 


0.27315E 


01 


0.27316E 


01 


0.27316E 


01 


0.27316E 


01 


0.27316E 


01 


0.27316E 


01 


0.27316E 


01 


0.27316E 


01 


0.27316E 


01 


0.27316E 


01 


0.27316E 


01 


0e27317F 


CI 


0.27317E 


01 


0.27317E 


01 


0.27317E 


01 


0.27317E 


01 


0.27317E 


01 


0.27317E 


01 


0.27317E 


01 


0.27317E 


01 


0.27 31 7E 


01 


0.27317E 


01 


0.27317E 


01 


0.27317E 


01 


0.273i7E 


01 


0.27317E 


01 


0.27317E 


01 


0.27317E 


01 


0.27318E 


01 


0.273L8E 


01 


0.27318E 


01 


0.27318E 


01 


0.27318E 


01 


0.27318E 


01 


0.27318E 


01 


0.27318E 


01 


0.27317E 


01 


0.27316E 


01 


0.27314E 


01 


0.27308E 


01 


0.27292E 


01 


0.27254E 


01 


0.27159E 


01 


0.26927E 


01 


0.26365E 


01 


0.25052E 


01 


C.22230E 


01 


0.17205E 


01 


0.12074E 


01 


C. 10389E 


01 


0. 10608E 


01 


0.18423E 


CI 
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PLUT OF GKUUP I / GKLUP 3 



0.163E 01 



7-AXlS 
0.206E 01 



0.2i^JE 01 



0.293E 01 



'.>.3J7F CI 



•J.C 

0.£>95b 01 

o.btn 01 

U.117E 02 
O.U6n 02 
0.17^>= 02 
0.2J51: 02 
0*23<»e: 02 
02 

0.292E 02 
3,32<:E 02 
a.35U 02 
v.idOF 02 
0.409E 02 
0.^3dE 02 
0.^6aE 02 
0.^^7l 02 
0.526= 02 
0.55^>E j2 
O.bdSb 02 
u.M^c 02 
0.643d 02 
0.t72E 02 
0.702fi 02 
0.731E 02 
0.7OJE 02 
U.TtJSE 02 
C.bl9r 02 
O.B'irdE 02 
O.B77= 02 
0,895t 02 
0.913E 02 
0.931*= 02 
0.9:>0F 02 
0.S66e 02 
O.Sdb2 02 
O.lOOt 03 
^1-^22 03 
J.lJ^P 03 
•J.lJfeL 0 3 
c.lCBt: Ui 

O^lO^c 03 
0.1 HE 0 3 
0.113c: 03 
0.115= 03 
O.ll/E 03 
0.119^ 03 
0.120r 03 
0.122k 03 
0.124t 0 3 
0.l2bF 03 
0.l27f: Oj 

oa2e6 03 

O.liOe 03 
0.l3lc 05 
0.132r 03 
0.13'tC 03 
0.135E 03 
J,ii7E 03 
0.138!: 03 
0.1396 Ci 
O.l^lE 03 
0.142=. 03 
0.144c C3 
0.l45c 03 
0. L4oS 03 
0.1486 03 
0.1496 03 
0.151E 03 
0.132E 03 
0.ib*»c 03 
0.l5S>f: 03 
G.l57e 03 
\i.l59c 03 
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PLOT OF GROUP 2 / GROUP i 



0.138E 01 



Z-AXIS 



0.20t)E: 01 



0.239F 01 



O.C 

J.5>tt3^ 01 
4.877E 01 
0«117£ OZ 
O.U6r 02 
0.175E J2 
0.2056 02 
0.2i4e 02 
0.263E 02 
0.2S^E 02 
02 

0.35lc 02 
O.SaOE 02 
0«A09b* 02 
0.4385 02 
0.46eE 02 
0.497B 02 
0.526H 02 
0,5:>5F 02 
Q.^83E 02 
G.614E 02 
0.643c 02 
0.672!: 02 
0.7025 02 
0.73i£ 02 
0.7605 02 
0.7Q9L 02 
J.B1S6 02 
0.e48C 02 
0.6775 01* 
0.«45t 02 
O.SlJc 02 
0.93iE 02 
0.9bOF 02 
J.S68S 02 
0.966E 02 
C.IOOE 03 
0.1D2E 03 
0.104E OJ 
0.106E 03 
O.IO&E 03 
0.109E 03 
O.lliE 03 
0.113F 03 
0.11?E 03 
0.117E 03 
0.119E 03 
0.120E ui 
0.122E 03 
0.124E 01 
0.125E 0 3 
0.127E 03 
0.128E 03 
0.130E 03 
0.131E 03 
0.132E Oi 
0.134E 03 
0.135E 03 
0.137E 03 
0.i38E 03 
^0.l39e 03 
'0«141E 03 
0.142E 03 
0.144E 03 
0.145E 03 
0.146E 03 
0.148E 03 
0.149E 03 
0.151E 03 
0.i52E 03 
0.154E 03 
0.155E 03 
0.157E 03 
0.159E 03 
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GROUP I / GROUP 2 



0.12327E 


01 


0.I2327E 


01 


0*12327£ 


01 


0.I2327E 


01 


0.12328c 


01 


0.I2328E 


01 


0.12329E 


01 


0.12329E 


01 


0.12329E 


01 


0.1233CE 


01 


0.12330E 


01 


0.12J31E 


CI 


0.12331E 


01 


0.12331E 


01 


0.12332E 


01 


0.12332E 


01 


0.12333E 


01 


0.1233^E 


CI 


C.12333F 


01 


0.12333E 


01 


0.1Z333E 


01 


0.12334E 


01 


0.l233^c 


01 


0.12334E 


01 


0.12335E 


01 


0.1233iE 


01 


0.12347E 


01 


0.12351E 


01 


0.1l8!)8t 


01 


O.aAAbE 


01 
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C.1.^327E 01 
C.1232flE 01 
0.12328E 01 
0.12329E 01 
C.12330E 01 
0.12331E 01 
0.12332E 01 
C.1233:E 01 
0.I2333E 01 
C.12333E 01 
0.12334E 01 
C.l233^^£ 01 
0.12337E 01 
0.12349E 01 
0.11A85E 01 



0.12327E 01 
C.i2328E CI 
0.12328E 01 
0.12329E 01 
0.i:33CE 01 
0.12331E 01 
0.12332E 01 
0.12332E 01 
0.12333E 01 
0.12333E 01 
C.12334E 01 
0.1233' E 01 
0.12339E 0\ 
0.12324E 01 
C.i2348E 01 



C.12327E 01 
0.12328E 01 
0.12328E 01 
C.12329E 01 
0.12330E 01 
0.12331E 01 
0.12332E 01 
0.12332E 01 
0.12333E 01 
0.12333E 01 
0. 1233<tE CI 
C.12335E CI 
0.12343E 01 
0.12217E 01 
0.15979E CI 



PLOT CF GROUP 1 t GROUP 



0.124E OL 



0.133E 01 



R-AXlS 

0.3 * 

o.<rgi!E 01 « 

O.^a^E 01 « 

0.d7/E 01 « 

Oi « 

U.U6E 0^ • 

u.ir^h 02 • 

0.2J6E Ml * 

O.Ji4E 0<! * 

0.26JE 0<J • 

0.292^ 02 * 

0.322E 02 • 

0.3i>lE J2 • 

o.3dOf 02 * 

0.<i09«! 02 • 

J.<.3dE J2 * 

0.^6(1': J2 * 

0.<»y7r: 02 • 

U.l'26k u2 * 

0.t>55-^ 02 • 

C. t?d3€ o2 • 
0.6l^E 02 • 
o.t^rJL 02 • 
J.fc72r 02 

0.7j2t C2 

0.7JIP 02 

J • fCiCc o2 < 

0.7d^r 02 ' 

a.rtmt 02 ' 

u.rt^d? 02 

o.*<77t 02 ' 

•j.t«)t)P 02 ' 

J.(»l3e 02 

0.«j3lE 02 

O.fj^UE 02 < 

J.*#oHf^ 02 ' 

0.slBo<: 02 

j.iooe OJ 

D. 102E 03 

O.lO^c OJ ' 
J. 106*^ Oj ' 

0. iJ^c OJ ' 

^•ili>: U3 ' 

j.ll*>E OJ 

0*U7n 0 3 

J. 11^5 Oj 

U.120E G3 

U.122F 03 

0.124E 03 

0.12;>^ OJ 

0.127E 03 

0.12dE 03 

O.liOE 03 

O.lilE 03 

0.132E 03 

0*134E 03 

0.135E Gi 

\>.137E OJ 

0.13BE J3 

0.lJ9t 03 

0.1<»IE 03 

0.1^2E 03 

J.1<»4E 03 

0.145E 03 

0.146E 03 

0.l4dE 03 

0.l49fi 03 • 

O.lblE 03 • 

0.152E 03 • 

03 • 

0.155E 03 ♦ 

0.1b7E 03 

0.159E 03 



POMCti OUTPUT IN KEGICnS 1*2*3 

0.13^U6E 10 0.77tt63E 09 0.32180E 09 



PQMEK density regions 1*2*3 (KU/tlTFR) 

0«lbl71t 03 J*133^aE 01 0.54903E 00 




JTED CVERRfcLAXATlCN FACTOR* 1.900 



